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CHAPTER  I 


INTRODUCTION 

Statement  of  the  Problem 

The  P-1 5 aircraft,  currently  entering  the  Air 
Force  inventory,  is  equipped  with  Pratt  and  Whitney  Air- 
craft F-100-PW-100  augmented  turbofan  engines  (8:2). 

The  F-100  engine  differs  from  older  engines  in  the  Air 
Force  inventory  in  two  major  respects:  1.  the  design 

incorporates  state-of-the-art  technology  and  2.  the 
engine  employs  modular  construction.  An  optimal,  field 
level  replacement  policy  for  the  engine  modules  is  needed 
(9). 

The  F-100  engine  is  composed  of  five  modules 
(11:48-49J.  As  shown  in  Figure  1.1,  they  are  the  fan 
module,  core  module,  fan  drive  turbine  module,  augmen- 
tor/exhaust  module,  and  the  gearbox  module.  In  addition 
to  the  modules,  a number  of  external  accessories  (e.g. , 
plumbing;,  wiring,  actuators,  probes,  valves,  pumps,  etc.) 
are  required  to  complete  the  engine  (1S:9-2A — 9-5).  When 
performing  maintenance  on  the  F-100  engine,  the  technician 
has  certain  sequences  he  must  follow.  These  sequences 
result  from  the  physical  order  of  the  modules  on  the 
engine.  For  instance,  to  remove  the  fan  drive-  turbine 


Figure  1.1 

F-100  Engine  Modules 

(Courtesy  Pratt  and  Whitney  Aircraft  Company) 
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module,  the  augmentor/exhaust  module  must  be  removed  first 
(18:3-39)-  Similarly,  in  order  to  remove  and  replace  the 
core  module,  both  the  fan  module  and  fan  drive  turbine 
modules  must  first  be  removed.  If  the  core  module  required 
replacement  and  there  were  nothing  wrong  with  the  fan  drive 
turbine  module,  the  latter  module  would  be  removed  solely 
to  facilitate  the  maintenance  on  the  core  module  and  then 
be  reinstalled.  There  may  be  opportunities,  however,  when 
it  is  advisable  to  replace  rather  than  reinstall  the  fan 
drive  turbine  module  when  it  is  already  removed  from  the 
engine.  To  look  at  the  extreme  case,  if  less  than  one 
operating  hour  remained  before  the  fan  drive  turbine  module 
were  due  mandatory  replacement,  re installation  of  the  same 
module  would  be  of  doubtful  merit..  At  the  opposite  extreme, 
if  several  thousand  hours  of  life  remained,  other  things 
being  equal,  replacement  of  the  module  would  be  unwise. 

Thus,  at  one  end  of  the  spectrum,  one  would  expect  to  incur 
excessive  maintenance-  manhours  and  downtime  because  the 
opportunity  to  accomplish  concurrent  maintenance  is  fore- 
gone; at  the  other,  much  useful  service  life  is  foregone. 
Somewhere  along  the  continuum  an  optimal  tradeoff  between 
maintenance  cost  and  lost  service  life  exists.  A policy 
which  optimizes  the  replace-not-replace  decision  for  the 
F-100  engine  modules  has  not  yet  been  developed  (7;9). 


H-W 
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Importance  of  the  Problem 

The  acquisition  cost  of  an  F-100  engine  is  approx- 
imately $1.7  million.  F-100  engine  replacement  parts  are 
correspondingly  expensive  (?)•  The  high  cost  of  the  F-100 
engine  and  its  sub-components  has  spurred  certain  efforts 
to  optimize  logistic  support  for  the  engine.  For  instance, 
the  Air  Force  Logistics  Command  (AFLC)  uses  a model  called 
MOD-METRIC  to  provide  improved  depot  management  of  inven- 
tory levels  for  both  the  basic  engine  and  its  modules 
(1 3:4-72). ^ MOD-METRIC  utilizes  historical  consumption 
data  and  forecast  flying  hours  in  order  to  compute  required 
engine  and  module  stocks.  The  end  objective  is  to  deter- 
mine the  smallest  stock  level  (and,  therefore,  cost)  which 
will  provide  required  support  (13:4-72).  MOD-METRIC  facili- 
tates decision  making  at  the  cf  V'  level.  Further 
reduction  in  cost  should  be  possible  if  the  base  level 
decision  on  when  to  replace  modules  can  be  facilitated  (9). 
Farther  development  of  the  example  given  previously  should 
clarify  this  point.  If  work  is  required  on  the  engine  core 
modulo  (which  requires  removal  of  the  fan  drive  turbine 


1 MOD-METRIC  is  an  extension  of  the  METRIC  (Multi- 
Koho] on-Tochnique-for-Recoverable-Item-Control ) model 
developed  for  the  USAF  by  the  RAND  Corporation  as  a method 
of  determining  stock  levels  for  recoverable  items  (1  '5:472). 
METRIC  addressed  a multi-echelon  (e.g.,  both  depot  and 
field  level),  multi-item  inventory  system.  MOD-METRIC 
extends  the  METRIC  model  to  include  an  indentured  inven- 
tory system.  The  F-100  engine  modules,  for  instance,  are 
an  indenture  of  the  basic  engine  (13:4-72). 
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module ),  under  what  situations  should  the  technician 
replace*  rather  than  simply  reinstall,  the  fan  drive 
turbine  module?  Under  the  current  operating  criteria, 
each  module  (except  the  augmentor/exhaust  module)  is 
given  a maximum  operating  time  (MOT)  at  the  end  of  which 
a time-change  is  required  (7).  In  our  example,  if  one 
hour  remained  on  the  fan  drive  turbine  module  until 
time-change  was  required,  the  technician  would  probably 
decide  to  change  the  fri  drive  turbine  module  now  rather 
than  reinstall  the  same  module  and  then,  one  operating 
hour  later,  remove  it  again.  These  factors  lead  to  the 
following  questions:  where,  in  the  continuum  of  pos- 

sibilities, is  the  bre.  joint  between  replace  and  not 
replace  which  will  result  in  minimal  cost  and  ?.  is  the 
effort  required  to  determine  the  breakpoint  worth  the 
savings  which  result  from  its  determination?  Determi- 
nation of  the  replace-not-replace  breakpoint  is  of 
considerable  current  interest  in  the  F-15  System  Program 
Office  (SPO),  at  Air  Force  Logistics  Command  (AFLC)  Head- 
quarters, at  Tactical  Air  Commc -ad  (TAC)  Headquarters,  and 
at  the  field  level  (7;9). 

Objectives 

The  objectives  of  this  thesis  were  to  develop  an 
algorithm  which  would  locate  the  economical 
replace-not-replace  breakpoint  for  the  five  modules  com- 
prising the  F-100-PW-100  engine  installed  in  the  F-^5 
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aircraft  and  to  investigate  sensitivity  of  this  breakpoint 
to  the  algorithm's  input  variables. 


Scope 

It  should  be  recognized  that  deciding  whether  or 
not  to  replace  an  unfailed  part  is  only  one  of  several 
similar  decisions  faced  by  the  maintenance  technician. 

For  example,  USAF  Technical  Order  2J-1-31  provides 
guidance  on  when  to  field  repair  and  when  to  return  to 
depot  components  which  have  already  failed  (31).  Tech- 
nical Order  2J-1-27  provides  guidance  on  when  to  minor 
overhaul  at  depot  level  and  when  to  major  overhaul  at  depot 
level  components  which  have  failed  (30).  This  thesis  will 
address  only  the  decision  on  whether  or  not  to  replace  an 
unfailed  moiule  at  the  field  level.  A statistical  model 
rather  than  inspection  criteria  will  be  used  to  anticipate 
failure. 


CHAPTER  II 


LITERATURE  REVIEW  AND  RESEARCH  QUESTIONS 

Air  Force  Jet  Engine  Maintenance  * 

The  Air  Force  engine  maintenance  program  consists 
of  two  echelons — Jet  Engine  Intermediate  Maintenance  (JEIM) 
and  depot  overhaul/repair  (2:vii).  JEIM  is  base  level 
repair.  Depot,  ovorhaul/repair  refers  to  extensive  tear- 
down  and  renewal  performed  at  either  an  Air  Logistics 
Center  or  contractor  facility. 

Periodic  maintenance  of  Jet  engines  consists  of 
scheduled  inspections  at  intermediate  level  and  scheduled 
major  overhaul  at  the  depot  level.  In  addition  to 
periodic  or  scheduled  maintenance,  unscheduled  maintenance 
can  and  does  occur.  Unscheduled  maintenance  requirements 
may  be  generated  for  either  the  base  level  or  depot  level. 
As  an  example,  catastrophic  failure  of  a component  would 
generate  an  unscheduled  maintenance  requirement  at  the 
base  and/or  depot  level.  Usually  unscheduled  maintenance 
is  more  costly  than  scheduled  maintenance  (10:225). 
Unscheduled  maintenance  is  basically  corrective,  whereas 
scheduled  maintenance  is  basically  preventive  (29:1-1). 


1 

By  catastrophic  failure  is  meant  failure  charac- 
terized by  sudden,  unexpected  damage  or  loss. 
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Standard  .let  engine  maintenance  cycle.  Figure  2.1  illus- 
trates the  standard  Jet  Sagine  Maintenance  Cycle.  The 
Standard  Jet  Engine  Maintenance  Cycle  is  the  most  common 
Ml;.  When  an  aircraft  engine  is  reported  operating  out- 
side of  established  performance  parameters  or  does  not 
meet  some  specified  inspection  criterion  an  unscheduled 
maintenance  action  is  generated.  A maintenance  technician* 
or  team  of  maintenance  technicians  is  dispatched  to  the 
aircraft  to  analyze  the  discrepancy.  Normally  the  engine 
will  be  analyzed  or  "trouble  shot"  to  determine  the  reason 
for  the  discrepancy.  Depending  upon  several  factors  (such 
as  estimated  repair  time,  type  of  failure,  inability  to 
determine  cause  for  failure,  availability  of  a repair 
asset  required,  etc.)  a decision  will  be  made  to  attempt 
repair  on  the  aircraft  or  remove  the  engine  from  the  air- 
craft and  take  it  into  the  shop.  Assuming  the  engine  was 

removed  and  is  operable,  the  next  step,  normally,  is  to 

-1 

perform  a test  cell  run  to  isolate  the  malfunction. 

Next,  the  engine  is  routed  through  the  jet  engine  repair 
shop  and  the  defective  components  replaced.  For  exces- 
sively damaged  engines  requiring  large  manhour  expenditures, 
a decision  may  be  made  to  ship  the  entire  engine  to  a depot 


1 

If  the  engine  is  not  operable,  the  malfunction 
is  known  with  a high  degree  of  certainty,  or  in  order  to 
save  time  the  engine  may  be  routed  directly  to  the  jet 
engine  repair  shop  (9). 
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Standard  lot  linginc  tttintcnancc  Cycle 
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Modular  Jet  l-ngine  Maintenance  Cycle 


Figure  2.1 

Typical  Jot  Engine  Maintenance  Cycle 
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overhaul  facility.  'Thus , the  engine  becomes  an  unscheduled 
requirement  at  the  depot. 

Engines  in  the  Air  Force  inventory  typically  have 
an  established  maximum  operating  time  (MOT).  On  reaching 
its  MOT,  an  engine  is  removed  from  the  aircraft  in  which 
it  is  installed  and  returned  to  the  depot  for  major  over- 
haul. Removal  on  reaching  MOT  is  a scheduled  maintenance 
action. 

F-100  engine  maintenance  cycle.  Figure  2.1  also  illus- 
trates the  F-100  engine  maintenance  cycle.  The  F-100 
engine  maintenance  cycle  differs  from  the  standard  jet 
engine  maintenance  cycle  because  of  the  modular  construc- 
tion of  the  F-100  engine  (11). 

Unscheduled  intermediate  level  maintenance  is 
essentially  the  same  until  the  F-100  engine  is  removed 
from  the  aircraft.  Normally,  a test  cell  run  is  made 

'l 

prior  to  engine  disassembly.  On  the  test  cell  run,  the 
defect  which  caused  engine  removal  is  isolated  to  a 
specific  module.  In  the  repair  shop,  the  defective  module 
is  removed  and  replaced  with  a serviceable  module  from  the 
stock  maintained  in  the  repair  shop  (1 3:4-73).  The  intent 
of  removing  and  replacing  modules  is  to  minimize  the  time 


1 

As  was  true  with  the  standard  jet  engine  main- 
tenance cycle,  the  engine  may  be  routed  directly  to  the 
repair  shop. 
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required  to  return  the  engine  to  serviceable  condition 
(7).  Engine  modularization  was  one  of  the  design  steps 
toward  minimizing  repair  cycle  time  on  the  P-15  aircraft 
(16).  Repair  of  moaules  is  scheduled  separately  from 
repair  of  the  entire  engine  and  does  not  delay  returning 
the  engine  to  serviceable  condition.  If  base  level  repair 
is  infeasible,  the  module  is  returned  to  the  depot  for 
overhaul  or  repair. 

Each  of  the  P-100  engine  modules  (except  the 
augment or /exhaust  module)  is  assigned  a maximum  operating 
time  (HOT).  When  a module  reaches  its  HOT,  the  engine  of 
which  the  module  is  a part  is  removed  from  the  aircraft 
and  brought  into  the  intermediate  level  repair  shop.  The 
module  which  reached  HOT  is  removed,  replaced  with  a 
serviceable-  module  from  the  chop  stock,  and  returned  to 
the  depot  for  overhaul. 

Optimal  Haintenance  Theory 

Jorgenson,  HcCall  and  Radner  (10:20-77)  in  a 
ItAIID  Corporation  report  entitled,  Optimal  Haintenance  of 
Stochastically  Failing  Equipment  provide  a comprehensive 
treatment  of  the  mathematical  determination  of  optimal 
maintenance  policies.  As  developed  in  their  report, 
maintenance  problems  may  be  divided  into  two 
classes — deterministic  and  stochastic  (10:1;.  Determi- 
nistic problems  are  those  where  the  requirements  and 
outcomes  of  every  maintenance  action  are  known  with 
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certainty.  For  stochastic  problems)  on  the  other  hand) 
the  requirements  and  outcomes  of  maintenance  actions  are 
random  in  nature  (10:1).  That  is,  the  amount  of  service 
life  produced  by  a unit  of  equipment  between  the  time 
when  a maintenance  action  is  performed  on  it  and  the  time 
of  failure  is  random  in  nature  rather  than  Irnown  before 
hand.  Jet  engine  failures  are  stochastic  in  nature;  in 
fact)  the  stochastic  nature  of  aircraft  engine  failures 
underlies  the  adoption  by  the  USA.F  of  the  actuarial  method 
of  predicting  aggregate  engine  failures  (29:1-1 — 1-4). 

The  question  formulated  in  the  statement  of  the 
problem)  earlier  in  this  thesis)  was  when  to  replace  an 
engine  module  which  had  already  been  removed  from  the 
engine  but  which)  itself)  had  not  yet  failed  or  reached 
MOT.  Replacement  and  not-replacement  were  the  only 
alternatives  considered.  This  limitation  is  reasonable 
since  the  field  level  maintenance  echelon  does  not  have 
the  option  of  repairing  an  engine  module  removed  by 
reason  of  accumulated  time$  but  must  return  it  to  the 
depot  (7;9). 

The  three  independent  variables  which  determine 
the  optimal  average  expenditure  per  unit  time  are 
replacement  cost  before  failure)  replacement  cost  after 
failure)  and  the  hazard  rate  (21:71-74).  Jorgenson, 

McCall  and  Radner  note  that: 

If  replacement  costs  more  after  a failure  than 
before)  in  the  absence  of  uncertainty  the  equipment 
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will  always  be  replaced  .just  before  it  fails.  For 
stochastically  failing  equipment , replacement  ,just 
before  failure  is  impossible  due  to  uncertainty 
about  when  failure  will  occur  [10:205  ] . 

For  the  stochastic  situation  it  may  or  may  not  pay  to 

replace  the  equipment  before  failure,  i.e.,  establish  a 

maximum  operating  time  or  opportunistic  replacement- 

policy.  Replacement  is  justifiable  if  two  conditions 

are  satisfied:  1.  the  time  to  failure  distribution  for 

the  piece  of  equipment  must  demonstrate  wearout,  and 

2.  it  must  cost  more  to  replace  the  piece  of  equipment 

after  failure  than  before  (10:221).  Replacement  prior  to 

failure,  however,  results  in  the  loss  of  some  unused  life. 

Thus,  the  optimal  decision  depends  on  the  tradeoff  between 

the  value  of  unused  life  and  the  cost  of  the  avoided 

failure  (10:207). 

A generally  used  technique  for  determining  if  an 
equipment  item  demonstrates  wearout.  is  examination  of  the 
hazard  rate  of  the  equipment  item  (24; 170-I71 Hazard 
rate  is  defined  as  the  ratio  of  the  number  of  failures 
occurring  in  a time  interval  to  the  number  of  equipment 
units  which  survived  until  the  beginning  of  the  interval, 
divided  by  the  length  of  the  time  interval  (24:161). 

Figure  2.2  illustrates  a typical  hazard  rate  function 
known  as  a "bathtub"  curve.  In  this  general  case  rlv- 
hazard  rate  initially  decreases  with  age,  remains  constant 
lor  a period  of  time,  and  th'-n  increases.  Th'-  equipment 
i tem  demonstrates  wearout  when  the  hazard  rate-  increases 
(24:171)  . 


DATA  HAZARD  RATE 


Figure  ?.? 

Typ  ical  Component  Failure  History 
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Even  if  the  equipment  item  exhibits  wearout, 
replacement  before  failure  would  still  not  be  justified 
unless  it  costs  more  to  replace  the  item  after  failure 
than  before  (10:223).  Replacement  cost*  either  before  or 
after  failure  would  generally  include  the  opportunity 
cost  of  the  downtime  while  the  equipment  is  out  of 
service  and  the  actual  dollar  replacement  cost  (10:224). 
Opportunity  cost  normally  is  a measure  of  lost  revenue 
(22:472-473).  National  defense,  the  product  of  military 
operations,  is  by  nature  not  priceable  (29:36-37).  Thus, 
ar  amount  of  lost  revenue  cannot  be  determined.  This  is 
not  to  say,  however,  that  the  opportunity  cost  of  incurred 
downtime  must  necessarily  be  ignored.  Jorgenson,  McCall 
and  Radnor  have  suggested  that  the  purchase  price  of  an 
equipment  item  be  amortized  over  its  total  expected  life. 

In  this  manner,  a cost  per  unit  of  operation  (e.g. , hour, 
cycle,  etc.)  can  be  determined  (10:224).  The  value  of 
the  amount  of  equipment  operating  time  foregone  by  some 
replacement  action  is  an  opportunity  cost. 

Intuitively,  one  would  expect  that  a replacement 
before  failure  would  result  in  less  downtime  than  a 
replacement  after  failure  (10: 222; . For  instance,  a 
replacement  after  failure  Is  generally  an  unscheduled 
maintenance  action.  As  an  unscheduled  action  is  unexpected, 
it  .generally  incurs  greater  queue  (waiting)  time  before 
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the  replacement  action  is  begun.  Jorgenson*  McCall  and 
Radner  (10:223)  point  out  that: 

...  a second  factor  determining  replacement 
costs  is  the  amount  of  resources  needed  to  perform 
the  action.  Usually,  more  resources  are  needed 
for  an  unscheduled  replacement  for  three  reasons: 

(1)  generally  a more  complex  maintenance  operation 
is  involved*  and  a :\iled  part  may  have  consider- 
ably lower  trade-in  value  than  an  unfailed  one; 

(?)  additional  resources  are  needed  to  repair 
parts  damaged  by  the  inservice  failure;  and 
fj)  it  is  often  necessary  to  transport  maintenance 
resources  to  the  failed  system. 

Under  certain  circumstances*  replacement  before 
failure  may  cost  less  for  at  least  one  additional  reason. 
To  use  the  F-100  engine  as  an  example,  when  core  module 
failure  or  removal  for  maximum  operating  time  (MOT)  is 
the  only  reason  for  engine  disassembly,  all  costs  asso- 
ciated with  removing  and  replacing  the  engine  in  the 
aircraft,  transporting  the  engine  to  and  from  the  repair 
shop,  and  removing  and  replacing  the  inlet  fan  module 
( in  order  to  gain  access  to  the  core  module)  are  sunk 
costs  against  the  core  module.  In  this  event,  if  the 
inlet  fan  module  were  replaced  with  a new  module  rather 
than  simply  being  reinstalled,  then  the  additional  engine 
removal  and  replacement,  engine  transportation,  and  module 
removal  and  replacement  costs  would  be  avoided. 

Replacement  of  the  inlet  fan  module  (or  any  other  module) 


The  inlet  fan  module,  gearbox  module,  fan  drive 
turbine  module  and  augmentor  module  must  be  remover,  to 
gain  complete  access  to  the  core  module  (18:3-2 — 3-7 1). 


given  that  its  removal  were  required  for  another  reason 
is  an  opportunistic  replacement  (10:224). 

In  summary*  the  following  factors  can  be  expected 
to  result  in  different  costs  for  replacement  before 
failure  and  replacement  after  failure: 

1.  Differences  in  downtime  required  for  the  replace- 
ment action. 

2.  Differences  in  resources  required  for  the  replace- 
ment action. 

3.  Opportunistic  replacement  when  removal  is  required 
to  support  mandatory  maintenance  of  another  component. 

Interviews  with  representatives  of  the  Product 
Support  Division*  Pratt  and  Whitney  Aircraft  Company  (27), 
the  prime  Air  Logistics  Center  (ALC)  for  the  F-100  engine 
(26),  and  personnel  of  the  Edwards  Air  Force  Base  F-15 
Joint  Test  Force  (3; 5; 12; 20)  were  conducted  to  determine 
if  it  is  possible  to  currently  measure  the  difference 
in  resources  or  downtime  for  replacement  before  and  after 
failure.  Although  the  personnel  interviewed  acknowledged 
chat  there  are  probably  such  differences,  the  data  to 
determine  such  differences  are  not  currently  available. 

The  difficulties  in  determining  resource  differences  can 
be  illustrated  by  examining  depot  overhaul  costs.  Whereas 
depot  overhaul  of  older  technology  engines  involves 
literally  complete  teardown  of  the  engine,  this  is  not 
true  for  the  F-100  engine.  Currently,  the  depot  overhaul 
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approach  for  the  F-100  engine  is  to  look  at  each  engine 
individually  and  only  replace:  1.  parts  which  have  met 

or  nearly  met  cycle  limitations  and  2.  failed  parts. 1 
Thus,  the  resources  required  to  overhaul  an  F-100  engine 
at  depot  are  dependent  not  only  on  the  reason  precipi- 
tating overhaul  (e.g. , MOT,  damage  beyond  field  level 
capability  to  repair,  or  opportunistic  replacement)  but 
also  on  the  age,  condition,  and  accrued  cycles  of  the 
engine  components.  .'Secondly,  data  available  on  F-100 
engine  depot  overhauls  are  very  limited.  The  San  Antonio 
ALC  (the  prime  ALC  for  the  F-100  engine)  received  the 
first  F-100  engine  for  overhaul  in  January  1975-  As  a 
result  of  the  uncertainty  over  what  maintenance  actions 
a typical  engine  will  require  and  the  lack  of  a developed 
data  base,  a distinction  cannot  as  yet  be  made  between 
the  cost  to  depot  overhaul  a failed  engine  or  module  and 
one  which  is  unfailod  (26). 

A similar  difficulty  exists  when  one  attempts  to 
determine  differences  in  resources  or  downtime  required 
for  intermediate  level  engine  maintenance.  Manhour  and 


1 

Cycle  in  this  context  is  used  in  a different 
sense  than  when  speaking  of  a maintenance  cycle  (see 
page  8 ).  Many  of  the  components  on  the  F-100  engine 
have  cycle  limitations,  where  a cycle  is  generally  con- 
sidered an  exercise  of  the  engine  throttle  from  idle  to 
an  advanced  power  setting  and  back  to  idle.  The  concept 
of  cycle  limitations  is  undergoing  much  discussion  at  the 
current  time  (7;9)  and  is  not  addressed  in  this  thesis. 
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clockhour  data  on  F-tOO  intermediate  level  engine  main- 
tenance are  obtainable  from  two  sources:  1.  the  Pratt 

and  Whitney  Aircraft  Company  F-100-FW-100  Qualitative, 
Quantitative,  Personnel  Requirements  Information  (QQPRI) 
(1?)  and  2.  from  data  maintained  by  the  Human  Factors 
Test  Office  of  the  F-15  Joint  Test  Force,  Edwards  Air 
Force  Base  (12).  Neither  the  QQPRI  nor  the  Human  Factors 
Test  Office  currently  makes  a distinction  between  the 
manhours  or  clockhours  to  perform  maintenance  on  unfailed 
modules  and  similar  requirements  for  failed  modules. 
First,  there  is  no  effort  at  this  time  to  make  such  a 
distinction  and  secondly,  it  is  questionable  if  such  a 
distinction  could  be  made  at  this  time  considering  the 
formative  state  of  available  data. 

Of  the  three  factors  which  can  result  in  a cost 
difference  between  replacement  of  a failed  module  and 
unfailed  module,  only  the  savings  through  opportunistic 
replacement  of  a module  which  is  already  removed  is 
tractable  at  the  current  time.  It  is  this  savings 
resulting  from  opportunistic  replacement  of  an  unfailed 
removed  module  which  is  explored  in  this  thesis.  As 
possible  savings  resulting  from  less  resources  or  less 
downtime  required  for  replacement  of  an  unfailed  module 
are  not  considered,  we  believe  that  total  savings  from 
opportunistic  replacement  are  understated.  The  effect 
of  uncertainty  about  the  total  savings  resulting  fro.m 
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opportunistic  replacement  was  explored  through  sensitivity 
analysis  and  is  discussed  later. 

Jorgenson*  lijCall  and  Radner  have  shown  how  to 
find  an  optimal  opportunistic  maintenance  policy  for  a 
system  composed  of  two  components*  one  of  which  is 
constant  hazard  (10:244-251).  The  Jorgenson,  McCall  and 
Radner  model  is  developed  in  terms  of  decision  rules. 

Over  the  interval  0<n£N  where  N is  the  maximum  oper- 
ating time  (MOT.)  and  n is  the  module  age  beyond  which 
opportunistic  replacement  is  worthwhile,  a module  is 
replaced  at  failure  in  the  interval  0<n,  replaced  at 
failure  or  opportunistically  in  the  interval  n£N,  and 
mandat orily  on  reaching  age  N.  The  values  of  n and  N 
which  will  result  in  least  cost  can  be  determined 
analytically  for  any  given  combination  of  cost  to 
replace  before  failure,  cost  to  replace  after  failure  and 
hazard  rate. 

Once  values  for  n and  N are  established,  they, 
in  effect,  form  a replacement  policy.  The  object  of 
this  thesis  is  to  develop  an  algorithm  which  will  enable 
managers  to  find  the  optimal  replacement  policy  for  the 
modules  of  the  F-100  engine.  For  the  F-100  engine, 
maximum  module  operating  times  (N)  have  been  established 
(7).  The  breakpoint  (n)  between  replacement  at  failure 
and  opportunistic  replacement  has  not  been  established. 
Thus , the  algorithm  developed  in  this  thesis  solves  only 
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for  the  optimal  value  of  n,  given  a fixed  value  for  N. 

In  the  literature  reviewed  on  optimal  preventative  main- 
tenance policy  determination  (6 : 271-283; “10: 205-268; 

19 : 229-249 ;21 : 61 -67; 35: 267-280)  a general  algorithm  which 
can  be  directly  applied  to  the  F-100  engine  module  waa 
not  found. 

Research  Questions 

1.  What  algorithm  can  bo  developed  for  a five  module 
system  which  1 for  any  combination  of  module  operating 
hours  and  hazard  rates  * will  determine  the  opt  imal 
opportunistic  replacement  policy? 

2.  How  sensitive  is  the  optimal  opportunistic 
replacement  policy  to  uncertainty  about  the  underlying 
failure  distribution? 

5.  How  sensitive  is  the  optimal  opportunistic 
replacement  policy  to  cost  estimate  uncertainty? 

4.  What  is  the  magnitude  of  the  savings  which  can 
bo  realized  by  an  opportunistic  replacement  poLicy  when 
compared  with  a replace  at  failure  policy? 


CHAPTER  III 


METHODOLOGY 


Introduction 

This  chapter  describes  the  method  used  to  answer 
the  research  questions.  The  procedures  used  to  determine 
hazard  rates  and  cost  data  are  described.  Next,  the 
structure  of  an  optimal  replacement  policy  for  the  five 
module  system  is  discussed.  Finally,  procedures  to 
determine  sensitivity  of  the  optimal  policy  to  failure 
rate  uncertainty  and  cost  uncertainty  are  described. 

Hazard  Rate  Models 

In  Chapter  II,  the  hazard  rate  was  described  as  a 
tool  for  determining  if  a component  exhibited  wearout  and 
as  an  input  to  determining  the  optimal  replacement  policy. 
This  section  will  discuss  how  hazard  rates  may  be  modeled. 

Ghooman  (24:160-170)  has  shown  how  to  develop  a 
hazard  rate  function  .from  failure  data.  If  the  time  scale 
in  a hazard  rate  graph  Is  divided  into  intervals,  the  data 
hazard  rate  for  each  interval  may  be  calculated  as  the 
fraction  of  components  surviving  until  the  beginning  of 
the  interval,  but  failing  during  the  interval,  divided  by 
the  interval  length.  Algebraically  the  data  hazard  rate 
is  defined  as 
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for  t . < tit.  + At.  where  Z.(t)  is  the  data  hazard  rate* 

is  the  number  of  survivors  at  the  beginning  of 

period  tj»  nft.^  + A t.)  is  the  number  of  survivors  at  the 

end  of  the  period  t^,  A t.  is  the  width  of  the  period  in 

time  units*  and  t is  the  total  observation  time.  As  an 
example,  if  there  were  113  survivors  at  the  beginning  of 
the  period*  2 4 failures  during  the  period,  and  the  period 
were  1000  hours 

Zd(t)  , 1m  . .0002124  (5.2) 


Although  the  hazard  rates  for  each  interval  may 
simply  be  plotted  on  a histogram,  Shooman  ( 24:185)  points 
out  that  in  order  to  generalize  from  sample  data  to  the 
population  of  similar  components  it  is  essential  to  fit 
the  failure  data  with  a mathematical  model.  Of  the  wide 
range  of  models  available,  Shooman  (24;195)  suggests  that 
the  piecewise  linear,  exponential,  and  Weibull  models  are 
sufficiently  inclusive  that  virtually  all  hazard  rates 
may  be  described  by  them.  Figure  3.1  illustrates  each  of 
those  models.  Included  with  each  illustration  is  the 
general,  algebraic  form  of  the  model. 
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Pitting  failure  data  with  a mathematical  model 
is  a two  step  process:  1.  choosing  the  appropriate 

model  and  2.  determining  model  parameters  (24:457). 

As  Shooman  (24:457)  indicates:  "The  initial  procedure 

in  the  choice  of  a m J'  1 is  to  plot  the  histogram  for 
Z^(t)  from  the  failure  data."  A judgement  is  then  made 

from  inspection  of  the  Z^(t)  function.  If  inspection  is 

not  adequate  to  determine  the  appropriate  model,  more 
powerful  analytical  and  graphical  techniques  are  avail- 
able (24:458-462;.  For  instance,  if  the  underlying  time 
to  failure  distribution  is  exponential  (i.e.,  the  hazard 
rate  is  constant),  a plot  of  the  number  of  operating 
hours  at  which  a component  failed  vs.  the  natural  loga- 
rithm of  ^ -j  where  i is  the  sequence  number  of  the  ith 

failure  and  N is  the  number  of  components  in  the  original 
population,  will  be  a straight  line  (24:459-460).  Similar 
techniques  are  applicable  to  the  Weibull  model  and  piece- 
wise linear  mode-1  ( 24 : 460-46  3 ; . 

bn  timet  ion  of  model  parameters,  once  the*  appro- 
priate model  has  been  selected,  can  be  accomplished 
through  either  least  squares  estimates,  moment  estimates, 
or  maximum-likelihood  estimates  (24:464).  Shooman  (24:464) 
recommends  the  maximum-likelihood  estimator  and  Jorgenson, 
McCall  and  Radner  (10:25^)  also  employ  the  maximum 
likelihood  estimator.  It  is  used  in  this  thesis  because 
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its  validity  and  applicability  is  widely  accepted  by 
researchers . 

Engine  Failure  Data 

Universe  description.  An  F-100  engine  prototype  was 
first  operated  in  flight  on  27  July  1972  (8:2).  Since 
that  time*  the  F-100  engine  and  its  modules  have  under- 
gone a number  of  modifications  to  resolve  various 
after-burner,  fuel  control,  and  other  component  malfunc- 
tions. As  a result,  the  F-100  engine  inventory  consists 
of  a number  of  different  configurations. 

Engine  Modules  Studied 

On  the  advice  of  the  Engine  Project  Office,  F-15 
System  Program  Office  (9),  the  population  of  F-100  engine 
modules  studied  was  limited  to  those  modules  originally 
installed  on  engines  serial  numbered  023  and  subsequent, 
with  the  exception  of  engine  number  050.  At  the  current 
time,  engines  serially  numbered  023  and  above,  except 
for  engine  number  050,  are  of  like  configuration  and  are 
the  operational  configuration. 

The  data  which  must  be  collected  in  order  to 
determine  a particular  module’s  (e.g. , fan  module,  core 
module,  etc.)  hazard  rate  are  the  times  at  failure  of 
those  modules  which  have  failed.  Time  at  failure  is  a 
discrete,  infinite,  ratio  level  random  variable.  Time 
at  failure  data  are  available  through  two  sources: 
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1 . the  Engine  Status  Reporting  System,  D024,  for  engines 
and  modules  assigned  to  operational  units  and  2.  through 
the  F-15  System  Program  Office  for  engines  and  modules 
assigned  to  Air  Force  and  Contractor  Research,  Develop- 
ment, Test  and  Evaluation  (RD'i'&E;  units  (9). 

Reporting  instructions  for  the  Engine  Status 
Reporting  System  are  contained  in  Air  Force  Manual  400-1, 
"Selective  Management  of  Propulsion  Units"  (54-:4-5)* 

The  source  document  utilized  in  recording  the  collection 
of  propulsion  unit  data  is  the  Air  Force  Form  1 554-  (54:5). 
Subsequent  to  completion  by  responsible  personnel  at  an 
Air  Force  Base,  AF  Forms  1554  are  key  punched  and  the  data 
transmitted  by  Automatic  Digital  Network  ('AUTODIN;  to 
Oklahoma  City  Air  Logistics  Center  (OCALC/ACDT; , Tinker 
Air  Force  Base,  Oklahoma  where  it  is  monitored,  processed 
and  maintained  ('54:12;.  Engine  or  module  operating  time 
at  failure  and  serial  number  are  specifically  collected 
by  this  system. 

The  Base  Engine  Manager,  designated  in  accordance 
with  Air  Force  Manual  400-1,  is  responsible  for  auditing 
and  controlling  AF  Forms  1554  submitted  from  his  base 
■'’54:  '0;.  He  develops  local  procedures  in  accordance 
with  which  an  initial  check  of  the-  accuracy  of  AF  Form 
'■  5 54  data  is  made-.  A second  check  on  data  validity  is 
accomplished  by  ^dit  routines  within  the  DO 24  system  ( 1 4 : i 9 ) - 
Finally,  at  the  end  of  each  month,  OCAMC/ACDT  provides 
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the  Base  Engine  Manager  with  a reconciliation  listing 
(D024AIH1A)  which  the  engine  manager  compares  with  his 
original  AF  Forms  "1 534-  and  verifies  for  accuracy  (34:5). 

As  discussed  above,  data  on  RDT&E  engines  are 
available  through  the  F-15  System  Program  Office.  The 
data  used  in  this  thesis  were  maintained  in  mechanized 
form  by  the  F-15  Join*;  Test  Force  at  Edwards  Air  Force 
Base,  California.  Although  the  data  are  not  subject;  to 
validity  checks  as  visible  as  those  provided  by  the  D024 
system,  F-15  System  Program  Engine  Project  Office  per- 
sonnel express  confidence  in  the  validity  of  their  data 
(9). 

Certain  engine  and  module  removals  are  precipi- 
tated by  events  other  than  engine  or  module  failure. 
Specifically,  it  was  necessary  to  purge  the  collected 
data  of  removals  to  facilitate  other  maintenance  on  an 
aircraft,  or  removal  caused  by  other  management 
decisions.  The  Engine  Status  Reporting  System  includes 
codes  to  distinguish  between  reasons  for  removal. 
Computer  programs  have  previously  been  developed  to  sort 
engine  failure  data  by  removal  code  (14:78-82).  For 
RDT&E  engines,  the  data  were  screened  utilizing  the  same 
programs  developed  for  D024  data  after  reformatting. 

As  the  F-100  engine  is  just  entering  the  inven- 
tory, total  operating  hours  and  the  total  number  of 
module-  removals  arc  quite  Low.  An  initial  computer  tape 
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containing  all  D024  reporting  trcuasacrions  for  F-100 
modules  through  30  March  1975  contained  only  one  reported 
module  failure.  A subsequent  tape  containing  data  through 
June  1975  revealed  three  additional  removals.  Data 
obtained  from  the  F-15  Joint  Test  Force  were  more  produc- 
tive though  still  containing  a limited  number  of  data 
points.  The  total  number  of  failures  and  maximum  time 
removals  by  module  type  were  as  follows: 

Table  3.1 

Maximum  Time  Removals  and  Failures  by  Module 

p 

Type  Reported  by  F-1 3 Joint  Test  Force 


Module  Type 

Failures 

Maximum  Time  Removals 

Fan 

16 

2 

Core 

5 

8 

Fan  Drive  Turbine 

12 

4 

Gearbox 

8 

6 

We  believe  the',  the  number  of  reported  module  failures  and 
maximum  operating  time  removals  ait  insufficient  to 
reliably  determine  the  underlying  hazard  rate.  Shooman 
f ^4: 4^7 ) , for  instance,  notes  that  statistical  techniques 


'The  augment or  module  is  not  assigned  a MOT  and 
for  this  reason  was  not  examined. 

p 

Methods  used  to  screen  the  raw  data  on  module 
transactions  are  described  in  Chapter  V. 
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of  analysis  begin  to  be  of  significant  benefit  when  at 
least  20  components  have  been  life  tested.  None  of  the 
modules  listed  in  Table  3.1  experienced  a total  number 
of  removals  (failures  plus  maximum  time  removals)  greater 
than  18.  However,  data  which  was  available  was  analyzed 
to  determine  a best  initial  estimate  of  the  hazard  rate 
function  for  each  module.  Hazard  rate  functions  thus 
determined  were  used  to  find  an  initial  optimal  value  of 
n (the  breakpoint  between  replace  at  failure  and  oppor- 
tunistic replacement)  and  to  investigate  sensitivity  of 
this  value  to  changes  in  the  hazard  rates. 

Cost  Factors  and  the  Structure  of  an  Optimal  Policy 

As  previously  pointed  out,  in  addition  to  knowing 
the  hazard  rate,  one  must  also  know  the  relative  costs 
of  replacement  before  failure  and  replacement  after 
failure  tr  determine  if  replacement  before  failure  is 
warranted.  In  Chapter  II,  the  factors  which  could  be 
expected  to  result  in  less  cost  for  replacement  before 
failure  were  found  to  be:  1.  fewer  required  resources, 

2.  less  downtime  for  replacement  and  5.  opportunistic 
replacement  when  an  equipment  item  was  already  removed. 
Further,  in  Chapter  II,  the  area  of  interest  for  this 
thesis  was  delimited  to  determining  an  optimal  policy 
for  opportunistic  replacement.  As  the  reader  will  recall, 
this  limitation  was  occasioned  by  the,  as  yet,  unsettled 
state  of  the  data  on  resource  costs  for  F-100  engine 
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maintenance  actions.  Opportunistic  replacement  policies 
have  been  explored  by  Jorgenson*  McCall  and  Radner 
(10:123-126)  and  others  (6 : 271 -283 ; 21 : 61-71 ; 38: 267-270), 
the  most  complete  treatment  being  by  Jorgenson,  McCall 
and  Radner.  A specific  approach  which  recognized  the 
savings  realizeable  through  replacement  of  a component 
which  has  already  been  removed  to  facilitate  other  main- 
tenance was  not  found  in  the  literature.  Such  an  approach 
is  developed  here  as  an  extension  of  the  Jorgenson,  McCall 
and  Radnor  model. 

Following  Jorgenson,  McCall  and  Radner  (10:244-251) 
one  may  divide  the  cycle  from  depot  overhaul  to  the  next 
subsequent  depot  overhaul  of  a single  module  into  two 
regions  as  illustrated  in  Figure  3.2. 


Replace  At  Failure  Only 


Replace  At  Failure 
Or  Opportunistically 


► 

0 


Region  A 


«+• 

n 

Module  Age 


Region  B 


N 


Figure  3.2 

Module  Depot  Overhaul  Cycle 
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In  region  A (i.e.,  age  less  than  n)  a module  would  be 
replaced  and  returned  to  the  depot  only  at  failure.  In 
region  B,  a module  would  be  replaced  both  at  failure  and 
opportunistically.  Opportunistic  replacement  is  meant 
to  be  the  replacement  of  a module  when  it  is  removed  to 
facilitate  maintenance  but  is  still  operable.  N is  the 
maximum  operating  time  where  modules  are  mandatorily 
removed  from  service.  For  modules  of  age  less  than  nt 
opportunistic  replacement  before  failure  costs  more  than 
reinstalling  the  same  module.  For  modules  of  age  greater 
than  n,  opportunistic  replacement  before  failure  costs 
less  than  reinstalling  the  same  module. 

The  mechanics  of  Figure  3.2  may  be  understood  by 
considering  the  underlying  costs  involved.  During  a 
single  cycle  from  depot  to  overhaul,  certain  costs  are 
always  incurred  and  others  may  or  may  not  be  incurred. 
Conceptually,  the  problem  is  similar  to  the  analysis  of 
fixed  and  variable  costs  discussed  in  economics 
''22:463 -46C,).  For  a single  overhaul  cycle,  the  fixed 
costs  are: 


1 

Specifically,  the  module  would  be  returned  to 
the  depot  only  if  it  experienced  a failure  beyond  the 
capability  of  an  intermediate  level  maintenance  facility 
to  repair.  When  gathering  data  on  module  times  at 
failure  in  support  of  this  thesis,  only  failures  precipi- 
tating a return  to  a depot  level  facility  were  counted. 


1.  The  cost  to  pack  a module  for  shipment  to  the 
depot  and  to  unpack  a module  upon  receipt  at  the  inter- 
mediate level. 

2.  Transportation  cost  from  the  intermediate  main- 
tenance level  to  the  depot  and  back. 

J.  The  cost  to  overhaul  the  module  at  depot.  As 
will  be  explained  in  Chapter  V this  cost  will  always 
be  assumed  to  be  the  same.  Further , this  cost  includes 
packing  and  unpacking  costs  at  the  depot. 

In  Chapter  II,  it  was  argued  that  how  one  assigns 
the  costs  to  remove  and  replace  the  engine  in  the  air- 
craft, transport  the  engine  to  and  from  the  intermediate 
maintenance  shop,  and  remove  and  replace  a module  on  the 
engine  depends  on  the  circumstances  causing  module 
removal.  Given  that,  for  instance,  the  fan  module  failed, 
then  these  costs  would  be  attributable  to  the  fan  module. 
On  the  other  hand,  if  the  core  module  had  failed  and 
removal  of  the  unfailed  fan  module  were  required  for  core 
module  repair,  then  these  costs  would  be  attributable  to 
the  core  module.  Subsequent  replacement  of  the  fan 
module,  given  that  the  engine  were  already  removed  from 
the  aircraft,  transport  e.'.  to  the  intermediate  maintenance 
shop,  and  disassembled,  would  not  generate  additional 
costs  to  perform  these  actions.  The  expected  engine 
removal  and  replacement,  transportation,  and  disassembly 
and  reassembly  costs  attributable  to  the  fan  module  depend 


on  the  probability  of  the  fan  module  being  replaced  either 
at  failure  or  alternately  upon  it  being  replaced  oppor- 
tunistically. 

One  may  note  by  reference  to  Figure  3 that  the 
probabilities  of  replacement  at  failure  and  opportunistic 
replacement  could  be  expected  to  bear  some  relationship 
to  the  sizes  of  region  A and  region  B.  Although  the 
specific  relationship  remains  to  be  examined*  it  does  not 
seem  unreasonable  that  as  region  A shrinks  in  relation  to 
region  B,  the  probability  of  opportunistic  replacement  of 
the  fan  module  would  increase.  Expected  expenditures  on 
engine  removal,  re installation,  transportation  and  disas- 

t 

sembly/assembly  attributable  to  the  fan  module  are  equal 
to  the  dollar  cost  to  perform  these  tasks  multiplied  by 
the  probability  of  incurring  them.  Thus,  as  region  B 
increases  in  proportion  to  region  A,  expected  expenditures 
on  engine  removal,  reinstallation,  transportation,  disas- 
sembly and  assembly  attributable  only  to  the  fan  modulo 
would  decrease. 

There  is,  however,  a penalty  for  increasing  the 
size  of  region  B.  In  region  A the  fan  module  is  replaced 
only  if  it  fails.  In  region  B the  fan  module  is  replaced 
for  failure  or  opportunistically.  As  region  B increases 
moves  towards  zero)  the  probability  of  early  oppor- 
tunistic replacement  increases  and  one  would  expect  that 
cycle  length,  i.e.,  the  mean  time  at  removal  (UTAH) 
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would  decrease.  Figure  3.3  illustrates  the  general  rela- 
tionship between  cost  and  cycle  length.  As  n moves  from 
N toward  0,  total  cost,  which  is  composed  of  fixed  cost 
and  variable  cost,  would  decrease.  The  cycle  length  over 
which  we  would  spread  these  costs,  however,  would  also 
decrease 1 

Thus,  it  is  not  immediately  obvious  if  a decrease 
in  total  cost  per  overhaul  cycle  is  advantageous  or  not. 
One  must  look  further  and  determine  cost  per  operating 
hour  (10:247).  The  optimal  value  of  n would  be  that  n 
which  results  in  minimum  cost  per  operating  hour. 

Following  our  discussion  thus  far,  suppose  that 
total  cost  is  some  function  f(n)  and  that  cycle  length 
is  also  a function  g(n);  both  f(n)  and  g(n)  are  monotonic 
and  increasing  in  the  interval  (0,  N>.  Let  h(n)  be  cost 
per  hour  as  a function  of  n where 


Following  accepted  methods  of  marginal  analysis  ( 18: 50-83 )» 
the  point  of  minimal  total  cost  per  operating  hour  could 
be  found  by  setting 

and  solving  for  the-  value  of  n which  would  make  this 
relationship  true.  Analytic  development  of  equation  (3.3) 
and  an  algorithm  for  solving  for  n are  discussed  in 


Cost/Hr.  Cycle  Length  Cost/Cycle 

(S  Per  Hr)  (Hr)  (I) 


Figure 

Cost  and  Cycle  Length  As  A Function  of  n 
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Chapter  IV,  Matheaatics  of  the  Algorithm. 

Sensitivity  of  the  Optimal  Policy  to  Unce^ainties  About 
Hazard  Rate  Function  I^arameters  and  Cost  Factors 

In  practice  the  parameters  of  the  hazard  rate 
function  are  not  generally  known  with  certainty  (10:?$y). 
furthermore,  the  total  cost  under  an  optimal  policy  may 
not  differ  significantly  from  the  cost  under  a 
replace-at-failure  policy.  As  an  example,  Jorgenson 
(10:228)  described  a system  whose  optimal  policy  produced 
only  a 2.5  percent  savings  over  the  replace-at-failure 
policy.  Therefore,  before  implementing  an  optimizing 
policy,  it  would  be  worthwhile  to  determine  the  policy’s 
sensitivity  to  uncertainty  in  hazard  rate  parameters  and 
cost  parameters. 

In  order  to  determine  sensitivity  of  the  optimal 
value  of  n and  the  optimal  policy  cost  to  uncertainty 
about  the  hazard  rate  parameters,  all  input  variables 
except  one  of  the  hazard  parameters  were  held  fixed  while 
that  parameter  was  varied  over  a range  on  either  side  of 
the  initial  estimate.  This  procedure  was  then  repeated 
for  each  of  the  other  parameters.  Sensitivity  to  cost 
uncertainty  was  also  examined  in  the  same  manner  by 
varying  one  cost  input  at  a time  while  all  other  inputs 
were  held  constant. 

finally,  a comparison  between  cost  of  operating 
the  system  until  failure  and  the  cost  of  optimal  replace- 
ment was  made  to  determine  the  significance  of  savings, 
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if  any. 

Sunnary  List  of  Assumptions 

1.  The  engine  failure  data  collected  for  this  thesis 
are  valid.  For  operational  engines  this  assumption  is 
based  on  an  examination  of  the  process  by  which  data  are 
generated,  recorded,  and  accumulated  as  described  earlier 
in  this  chapter.  For  RDT&E  engines,  this  assumption  is 
based  on  discussions  with  the  Engine  Project  Office,  F-15 
System  Program  Office. 

2.  The  hazard  function  of  each  of  the  five  modules 
is  independent  of  the  other  modules.  This  assumption  is 
based  on  discussion  of  the  two-component  system  by  Jorgenson, 
McCall  and  Radner  (10:244-25'!;. 

3.  Transportation  costs  and  depot  overhaul  costs  for 
modules  returned  to  the  overhaul  facility  are  assumed  to 
be  the  same  whether  the  module  has  failed  or  not. 

Summary  List  of  Limitations 

1.  The  algorithm  developed  in  this  thesis  to  deter- 
mine an  optimal  replacement  policy  is  applicable  only  to 
a five  component  system.  Similar  methodology  can  be  used 
in  other  systems.  Further,  the  only  decision  the  policy 
will  facilitate  is  whether  or  not  to  replace  a component 
which  has  not  yet  failed. 

2.  Sonora  1 izat ions  derived  from  the  data  gathered 
by  this  research  can  bo  made  only  to  the  population  of 
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F-100  engine  modules  originally  installed  on  F-"IOO  engines 
023  through  049  and  engine  051 » and  to  subsequent  modules 
meeting  the  same  design  specifications. 

3.  Opportunistic  policies  for  unfailed  modules  that 
are  not  removed  to  facilitate  other  maintenance  were  not 
considered. 

4.  Conceptually*  it  is  optimal  economically  to  solve 
for  n and  N simultaneously  but  in  this  thesis  N is  given. 
Logically,  In  practice  we  grope  toward  an  appropriate  N 
through  careful  observation  of  failures. 


CHAPTER  IV 


MATHHIATICS  OP  THE  ALGORITHM 


Introduction 

This  chapter  presents  the  mathematics  of  the 
algorithm  used  to  find  an  optimal  value  of  n for  a given 
set  of  hazard  function  parameters  and  cost  factors.  In 
essence,  the  algorithm  is  an  extension  of  a model  developed 
by  Jorgenson,  McCall  and  Radner  (10:?44— 251 ) for  a two 
component  system.  The  current  formulation  assumes  con- 
tinuous underlying  time  at  failure  distributions,  permits 
the  second  component  (in  this  case  the  core  module)  to  have 
a general  hazard  function,  and  solves  for  n;  the  Jorgenson, 
McCall  and  Radner  model  on  which  it  is  based  assumed 
discrete  distributions,  restricted  one  system  component 
to  a constant  hazard,  and  solved  simultaneously  for  n 
end  N. 

The  order  of  presentation  is  as  follows:  first, 

mathematical  expressions  are  derived  to  calculate  expected 
cost  per  cycle;  second,  expected  cycle  length  is  addressed; 
third,  calculation  of  conditional  probabilities  is  dis- 
cussed; fourth,  the  method  of  incorporating  core  module 
age  is  introduced  and  finally,  the  minimization  technique 
is  discussed. 
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Table  4.1  indicates  which  modules  must  be  removed 
to  facilitate  removal  of  other  modules.  It  can  be  seen 
that  only  core  module  and  fan  drive  turbine  module  removal 
precipitate  removal  of  other  modules.  It  will  be  recalled 
that  our  interest  is  in  when  to  replace  an  unfailed  module 
which  has  already  been  removed  to  facilitate  maintenance 
on  another  module.  It  can  be  seen  from  Table  4.1  that  the 
fan,  fan  drive  turbine,  augment or/exhaust  and  gearbox 
modules  under  certain  circumstances  will  require  prior 
removal  to  facilitate  maintenance  on  another  module.  The 
augment or/exhaust  module,  however,  does  not  have  an  estab- 
lished depot  overhaul  interval  (7).  All  repair  is 
accomplished  at  the  field  level  and,  from  the  standpoint  of 
the  algorithm  developed  in  this  chapter,  the 
augmentor/exhaust  module  will  not  be  considered.  The  fan 
module,  fan  drive  turbine  module  and  gearbox  modules  are 
removed  to  facilitate  maintenance  on  the  core  module.  Thus 
we  need  only  explore  when  to  replace  unfailed  fan  modules, 
fan  drive  turbine  modules  and  gearbox  modules  when  they 
are  removed  to  facilitate  maintenance  on  the  core  module. 

The  algorithm  is  developed  below  in  terms  of  the  fan  module 
oat  is  equally  applicable  to  the  fan  drive  turbine  and  gear- 
sox  modules. 
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Opportunistic  replacement  before  failure  is 
warranted  only  if  1.  opportunistic  replacement  costs 
less  than  reinstalling  the  module  and  2.  the  module  has 
an  increasing  hazard  rate  in  the  region  where  opportunistic 
policies  are  considered.  Whether  or  not  a module  satisfies 
the  second  criterion  is  determined  by  examining  its  hazard 
rate  as  discussed  in  Chapters  III  and  V.  The  algorithm 
developed  in  this  chapter  would  be  exercised  only  if  the 
hazard  criterion  were  met. 

Determination  of  Expected  Cost 

Figure  3.2  from  Chapter  III  is  redrawn  here  for 
ready  reference.  The  reader  will  recall  that  in  region 
A,  the  fan  module  is  replaced  and  returned  to  the  depot 
only  if  it  sustains  a failure  requiring  depot  overhaul. 

In  region  B,  the  fan  module  is  replaced  and  returned  to 
the  depot  if  it  either  fails  or  reaches  its  maximum  oper- 

'l 

ating  time  or  if  the  core  module  is  removed  for  any  reason. 


1 

It  should  be  noted  that  the  possibility  of  inter- 
mediate l«*vei  reparable  failures  exists  in  both  regions. 
Thus,  it  is  necessary  to  screen  data  from  which  hazard 
parameters  are  determined  for  the  fan  module,  fan  drive 
turbine  module  and  gearbox  module  in  such  a manner  that 
only  depot  reparable  failures  are  included.  The  simpli- 
fying assumption  is  made  here  that  all  failures  for  module 
age > n are  depot  reparable  only.  Scrictly,  this  is  not 
likely  to  be  true.  The  direction  of  probable  bias  resul- 
ting from  this  simplification  is  discussed  in  Chapter  VI. 
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Module  Depot  Overhaul  Cycle 
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N 

n 

Q(to) 

Q.(  fo) 

Ct(fc) 

q(f) 

qi(fc) 

qi(i'cj 

QFAN(t 


Age  of  ue  fan  module  at  which  mandatory 
replace  -lent  of  the  fan  module  occurs  given 
survival  until  that  age, 

i 

Earliest  age  of  the  fan  module  at  which 
opportunistic  replacement  is  permitted: 
n is  the  leading  edge  of  region  B, 

Probability  that  an  engine  requires  field 
level  replacement  of  the  fan  module  and  the 
core  module  simultaneously  between  n and  N 
given  survival  of  both  until  n, 

Probability  that  an  engine  requires  field 
Level  replacement  of  the  fan  module  before 
the  core  module  between  n and  N given 
survival  of  both  until  n, 

Probability  that  an  engine  requires  field 
level  replacement  of  the  core  module  before 
the  fan  module  between  n and  N given  survival 
of  both  until  n. 

Probability  of  fan  replacement  before  age  n, 

Probability  that  an  engine  requires  field 
level  replacement  of  the  fan  module  and 
core  module  simultaneously  in  the  i-th 
interval  between  n and  N given  survival  of 
both  until  n, 

Probability  that  an  engine  requires  field 
replacement  of  the  fan  module  first  in  the 
i-th  interval  between  n and  N given  survival 
of  both  until  n, 

Probability  that  an  engine  requires  field 
replacement  of  the  core  module  first  in  the 
i-th  interval  between  n and  N given  survival 
of  both  until  n, 

I ^ 

Pan  JFan 

Probability  of  fan  module  failure  between 
fan  module  age  and  fan  module  age  t^ 

given  survival  until  fan  module  age  t . , 

<J 
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QCORE(tk  It,  ) 
core 1 Jcore 


Probability  of  core  module  failure  between 
core  module  age  t..  and  core  module  age  t^ 

given  survival  until  core  module  age  t., 

tJ 

f w f t ) Pan  module  failure  probability  density 

ran  function, 

f Uj  Gore  module  failure  probability  density 
" ' function, 


PANREP 

FANDEP 

FANPACK 

FANSHIP 

FANMOT 

GOREMOT 


Fan  module  field  level  replacement  cost, 

Fan  module  depot  overhaul  cost, 

Cost  to  package  the  fan  module  for  shipping 
and  unpackage  upon  receipt, 

Round  trip  shipping  cost  for  the  fan  module 
between  the  field  level  repair  shop  and  the 
depot , 

Fan  module  maximum  operating  time, 

Core  module  maximum  operating  time. 


The  probability  of  fan  module  removal  for  failure  in  the 
interval  (0,  n)  is  simply 


<?(*) 


W*> dt 


(4.1) 


Consider  an  interval  of  width  6 t between  n and 
N where  at  « fei+l  “ fei  and  It  is  sufficiently  small  that 

only  one  of  the  possibilities  q^fc),  q^fcT),  or  qi(Fc) 

can  occur  in  the  interval.  During  the  i-th  interval, 
the  probability  that  the  fan  module,  the  core  module,  or 
both  are  removed  given  both  are  installed  at  the  same  time, 


both  are  the  same  age , and  given  survival  of  both  until 
n is 
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q^f  Uc|  *0  = q^fc)  + q^fc)  + q^fc).  (4.2) 

When  n<  L SN  - It,  then 

qi(fc;  = [l  - QPAN(ti|  n)]*  (4.5) 

[l  - QCORE(ti|  n)j« 

QFAN(ti+1|  tt)* 

QCORE(ti+1 | t.), 

and  when  N - I t < N and  N * COREMOT 

q^fc)  - [l  - QFAN(N  - I t|  n)]«  (4.4) 

[l  - QCORE(N  - I t|  n)]« 

QFAN(N|  N - I t) 

since  core  removal  would  then  be  certain.  Similarly,  if 
N - i t<t^N  when  N * FAJNMOT  then 

qt(fc)  - [l  - QFAN(N  - I t | n)]*  (4.5) 

[l  - QCORE(N  - I t | n)J* 

QCORE(N  | N - It) 

The  terms  q^(fc)  and  q^(fc)  are  developed  similarly.  When 
n < t . < N - It,  then 


1 

The  conditional  probabilities  developed  in  this 
chapter  may  not  be  intuitively  clear.  A more  complete 
mathematical  treatment  is  found  in  Appendix  H. 
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and, 

but 

Also 


but , 

but 

Then 


qi(fc)  = [l  - QFAH(tJ  n)]* 

[l  - QCOHE(ti+1 | n)]* 

when  N - ft  and  N = COREMOT  qi 

if  N = PANMOT  then 

r\L(rZ)  = [l  - QPAN(N  - ft  t | n)J. 

[l  - Q«OkE(N  | n)]  # 

, when  n < t^<  N - ft  t, 

q^Fc)  = [l  - QFAN(t.+,,  | n )j  • 

['  - QCORE(t.|  n)J. 
QCORE(ti+1|  t.) 

when  N - ft  t <t.  and  N = COREMOT 

q^fc;  = - QFAN(N  | n)J» 

[l  - QCOREfN  - ft  t | n) 
3i(fc ) = 0 if  N = PANMOT . 

, where  M - (N  - n)/  ft  t such  that  M is 

M 

Q(fc)  = £ 9j_(fc) 
i='1 

M 

Q(fcj  = E and 


(4.6) 


(fc)  = 0, 

(4.7 ) 


(4.8) 


(4.9) 

]• 

an  integer, 

(4.10) 


(4.  ) 
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Q(fc)  - L 

i=1 

It  should  be  noted  that  q^fc)*  q^Cfc),  and  q^fc)  are 

approximations  of  Q(fc),  ^jr  Q(fc),  and  ^ Q(fc) 

respectively.  The  accuracy  of  the  approximation  depends 
on  how  many  segments  the  interval  (n,  N)  is  broken 
into — that  is  on  the  size  of  * t.  A more  complete  dis- 
cussion of  approximation  accuracy  is  found  in  Appendix  F. 

FANPACK,  FANSHIP  and  FAJNDEP  were  previously  defined 
an  the  cost  to  package  and  unpackage  the  fan  module,  ship 
the  module  to  and  from  the  depot  and  accomplish  fan  module 
depot  repair.  FANPACK,  FANSHIP  and  FANDEP  are  experienced 
with  certainty  during  a single  overhaul  cycle.  FANREP, 
the  cost  of  field  level  removal,  however,  is  experienced 
only  if  the  fan  module  fails  in  the  interval  (0,  n)  or  is 
removed  before  the  core  module  in  the  interval  (n,  N). 

It  would  not  be  experienced  if  the  core  module  failed 
before  the  fan  module  in  the  interval  (n,  N).  Therefore, 
the  expected  cost  in  the  interval  (0,  N)  is 

K(cost  ) - (4.  Ij) 

FANPACK  + FANSHIP  + FANDEP 

+ FANREP  • jq(f ) + [l  - q(f)]*[Q(fc)  + 

-] 

We  will  use  the  convention  that  if  both  the  fan 
module  and  core  module  require  removal  at  the  same  time 
(e.g. , simultaneous  failure)  FANREP  is  chargeable  to  the 
fan  module. 


■2v»r*r 
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Determination  of  Expected  Cycle  Length 

The  expected  value  of  a discrete  probability 
distribution  is 

K(k)  = *k  = £ kKk)  (4.14) 

where  k is  a discrete  random  variable  and  P(k)  is  the 
probability  of  observing  the  value  k (57:148).  The 
equivalent  'repression  for  a continuous  distribution  is 


E(x)  = 


(4.15) 


where  x is  a continuous  random  variable,  f(x)  is  its 
density  function  and  a and  b are  the  limits  of  integration 
(57:195). 

In  the  interval  (0,  n),  the  probability  of  fan 
module  replacement  was 


<3  (f  ) 


rFan 


( t ) dt 


(4.1) 


and  in  the  interval  (n>  N)  the  probability  of  fan  module- 
replacement  is 

P( replacement  in  interval  (n,  N))  » (4.16) 


M 


[1  - q^f)]  •£[qi(fc)  + 9i(-fc)  + 9i(^] 
i=1 


'he-  expected  t ime-  at  re-placement  of  the  fan  module  then  is 


E(time  at  replacement)  « (4.1?) 

jF  * • fFao(t)  dt 

+ [ 1 - q(f )]• 

M-1  + t ) 

1?1  { — L'g~— - f)i(rc)  + Qi(fc)  ♦ q^Fc)]} 

+ [1  - q(f)]*[1  - QCORE(N  | 0)  ] S 

where  the  term 

[1  - q(f)]«[l  - QC0RE(N|  0)] 

i«  the  probability  of  fan  module  and  core  module  survival 
until  N.  The  expected  time  at  replacement  is  also  the 
expected  cycle  length  since  the  fan  module  is  returned 
to  the  depot  on  replacement. 


Conditional  Probabilities 

QFM(t.  t . ) and  QC0RE(tk  It.  ) are 

KFan  '■'Fan  *core  I ‘'core 

conditional  probabilities.  Using  a derivation  by  Papoulis 

(15:179),  the  conditional  probability  of  module  removal 

between  t^  and  t^  given  survival  until  t^  can  be  expressed 

as 

, F(tk)  - F(t.)  1Rs 

H removal  | survival  until  t^)  = — PCt'.';  ^ 


It  is  assumed  that  replacement  occurs  at  the  mid- 
point of  the  interval  (t^,  t^+/j). 
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where  F(t)  is  the  cumulative  distribution  function  evaluated 
at  t. 

Core  Module  Age  and  Difference  in  MOT's 

It  would  not  normally  be  true  that  both  the  fan 
module  and  core  module  would  have  the  same  number  of 
accumulated  hours.  To  illustrate  with  an  obvious  example, 
if  at  some  fan  module  age  t a new  (0  time;  core  module 
were  installed  on  the  engine,  the  difference  between  core 
module  and  fan  module  age  would  be  t hours.  As  core 
module  age  can  have  significant  effect  on  the  probability 
of  core  module  removal,  it  is  important  to  provide  for 
core  module  age  as  an  input  to  the  algorithm.  Further, 
although  at  the  current  time  both  the  core  module  and  the 
fan  module  have  the  same  maximum  operating  time  (riOT) 
this  will  not  necessarily  be  true  in  the  future.  A 
difference  in  fan  module  and  core  module  MOT  is  also  an 
important  input  t . the  algorithm.  To  illustrate,  if  the 
fan  module  and  core  module  were  of  the  same  age  at  the 
time  of  fan  module  opportunistic  replacement  and  the  core 
module  had  the  earlier  MOT.  the  fan  module  cycle  could  not 
last  longer  than  the  amount  of  time  remaining  on  the  core 
module.  Core  module  removal  at  its  MOT  given  survival 
to  its  MOT  is  a certainty.  Therefore,  we  would,  with 
certainty,  ship  the  fan  module  to  the  depot  on  achievement 
of  core  module  MOT  under  an  opportunistic  policy. 


l.JW  JlLim-TIFF  Z~' 


wrrK^mi^m 


IPW*1  .'•  1 " 
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The  following  method  will  he  used  to  incorporate 
core  module  age  and  a difference  between  fan  module  MOT 
and  core  module  MOT  into  the  algorithm.  Let 
FANTIME  B Age  of  fan  module 
CORETIME  s Age  of  core  module 
A time  5 FANTIME  - CORETIME 

and 

A MOT  B FANMOT  - COREMOT 

Then,  alternatively,  we  may  express  CORETIME  and  COREMOT 
in  terms  of  FANTIME  and  FANMOT  as 

CORETIME  = FANTIME  - A time  (4.19) 

and 

COHEMOT  = FANMOT  - A MOT.  (4.20) 

FANMOT  could  be  reached  first  if 

FANMOT  - FANTIME  < COHEMOT  - CORETIME 
which,  by  definition,  is  the  same  as 

FANMOT  - FANTIME  < FANMOT  - A MOT  - (FANTIME  - A time) 
or,  more  simply 

A time>  AMOT. 

Then,  if  A time  > A MOT 

N = FANMOT,  (4.21  ) 

and  if  A time  = A MOT 

N - FANMOT  = COHEMOT  (4.22) 

but  if  A time  < A MOT 


N = FANMOT  - A time  + A MOT 


(4.2J) 


5 4 


The  conditional  probability  QCOHE(t.  It.  ) is 

*core I “core 

read  as  the  probability  of  core  module  failure  between  core 
module  age  t^  and  core  module  age  t^  given  survival  until 

core  module  ago  t..  When  the  difference  between  PANTIME 

U 

and  COliETIMK  is  A time,  the  conditional  probability 

QCOIiE(t.  It.  ) can  be  expressed  in  terms  of  PANTIME 
Kcore ' “core 

as  QCORE(tv  - A time  I t.  - A time)  where  this  expres- 
*fan  1 Jfan 

sion  appears  in  Equations  (4.3)»  (4.4),  (4.5)»  C-4-* B ) * (4.7), 
(4.8),  (4.9)  and  (4.17). 

Minimization  Technique 

Expected  cost  per  hour  is  determined  as  the  ratio 
of  expected  cost  per  cycle  to  expected  cycle  length. 

Expected  cost  per  cycle  and  expected  cycle  length  are 
expressed  by  Equations  (4.13)  and  (4.17)  respectively. 
Minimization  of  this  ratio  is  tractable  through  numerical 
techniques  (10:123).  The  technique  used  in  this  thesis 
was  to  calculate  the  expected  cost  per  hour  as  the  value 
of  n was  stepped  in  10  hour  increments  from  0 to  250 
hours.  This  process  was  then  repeated  as  core  age  was 
increased  in  ten  hour  increments  from  zero  to  250  hours. 

For  each  value  of  core  age,  the  optimal  value  of  n was 
that  value  which  resulted  in  least,  expected  cost  per  hour. 
Ohoico  of  ten  hour  increments  for  n and  core  age  was  in 
large  measure  arbitrary.  Obviously,  one  could  make 
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the  interval  as  sma  as  desired.  Smaller  intervals, 
however,  exact  a penalty  in  computer  processing  time. 
Using  ten  hour  intervals,  1.?  processing  hours  were 
required  to  run  the  program  listed  in  Appendix  F in  the 
batch  mode  on  the  computer  available  at  Wright-Pat terson 
Air  Force  Base. 


CHAPTER  V 


DATA  ACQUISITION  AND  ANALYSIS 


Cost  Data 

Several  elements  make  up  the  costs  involved  in 
modular  engine  repair  actions  (5).  The  costs  addressed 
in  this  research  effort  were:  1.  module  depot  overhaul 

costs , 2.  base-to-depot  and  return  transportation  costs, 
3.  module  packing  and  unpacking  costs  and  4.  inter- 
mediate (field)  level  repair  costs.  Each  of  these  cost 
elements  will  be  discussed  separately  in  this  section. 
Detailed  calculations  and  the  data  used  are  found  in 
Appendix  B.  For  consistency,  all  base  level  repair  costs 
are  assumed  to  take  place  at  Edwards  Air  Force  Base, 
California. 

Depot  costs.  Overhaul  at  depot  level  is  accomplished  to 
return  a module  to  like-new  configuration.  Several 
factors  such  as  depot  material  cost,  depot  manhour  cost 
and  depot  overhead  are  elements  of  the  single  cost 
referenced  in  this  thesis  as  depot  overhaul  cost.  In 
Chapter  II  it  was  noted  that  depot  overhaul  cost  per 
engine  or  module  can  be  highly  variable.  Thus,  reduction 
of  depot  overhaul  cost  to  a single  figure  ignores  the 
considerable  uncertainty  inherent  in  this  cost;. 

36 


More 
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detailed  cost  data*  however*  is  not  yet  available.  For 
this  "eason,  the  officially  accepted  module  depot  over- 
haul costs  were  utilised  as  a starting  point  to  determine 
an  optimal  policy  and  to  investigate  sensitivity  of  the 
policy.  The  officially  accepted  module  overhaul  costs 
were  taken  from  San  Antonio  ALC/MAWWF  letter  of  12  March 
1975  entitled  Module  Depot  Overhaul  Costs  (7). 

Base  to  depot  transportation  costs.  Engine  modules  for 
the  F-100  engine  are  shipped  by  truck  and  classified  as 
machinery,  parts,  steel  f j.  Charges  for  this  movement 
are  based  solely  upon  weight.  AFLC/MMP  provided  a '''hurt 
showing  the  weight  of  each  module  when  properly  packed 
in  its  shipping  container  (7;.  A shipping  rate  schedule 
( Edwards  Air  Force  Base,  California  to  San  Antonio  Air 
Logistics  Command,  Texas;  which  showed  cost  per 
hundred-weight  was  supplied  by  the  Office  of  the  Chief  of 
Iransportation/DSTRP,  Kelly  Air  Force  Base,  Texas.  Multi- 
plication of  the  appropriate  rate  schedule  by  the  module 
shipping  vr.ight  gave  the  cost  of  one  way  movement  of  each 
individual  module-.  Since  overhaul  action  requires  a now 
module  to  be  sent  to  the  base,  as  well  as  shipment  of  the 
oLd  one  to  depot,  complete  transportation  cost  involves 
two-way  shipment  which  doubles  the  cost  thus  calculated. 

Intermediate  c fields  level  costs.  Several  costs  are 
incurred  at  field  level,  consisting  of  manpower  charges 
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and,  if  appropriate,  material  charges.  All  manpower  costs 
were  figured  based  upon  standard  team  sizes  of  four,  three 
or  two  technicians  as  required  for  the  action.  Costs  per 
clockhour  were  calculated  for  each  size  standard  team 
based  upon  wage  rates  found  in  AFM  1?7-101.  Three  sources 
were  used  to  obtain  estimates  for  the  clockhours  required 
for  each  task.  The  three  potential  sources  were:  1.  Quali- 

tative and  Quantitative  Personnel  Requirements  Information 
( QQPRI ) M'/j,  a document  provided  by  Pratt;  and  Whitney  to 
predict  manhour  requirements,  2.  time  and  motion  studies 
from  videotape  films  produced  at  Edwards  Air  Force  Base 
during  technical  order  validation  work  and  3.  subjective 
estimates  based  upon  a supervisor's  actual  experience  on 
the  F-1Q0  engine.  A comparison  of  the  three  estimates  is 
presented  in  Table  E.  "11.  The  QQPRI  figures  were  predic- 
tions made  based  upon  engineering  design.  Validation  of 
the  figures  is  as  yet  incomplete  but  preliminary  results 
wore  inconclusive  M Time  and  motion  studies  from 
videotape  films  were  found  to  be  poor  est imatos  since 
technical  order  validation  action  requires  considerable 
time  to  stop  and  document  difficulties  found.  Subjective 
estimates  based  upon  experience  were  declared  to  be  the 
mos ' consistent,  valid  figures  available  at  this  time 
'"1?;20;.  Since  these  estimates  reflect  100  percent 
utilization  of  personnel,  an  allowance  factor  of  1.67 
provided  by  the  F-15  JTF  at  Edwards  Air  Force  Base  was 
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applied  to  convert  these  estimates  to  expected  average 
performance.  Manhour  estimates  are  available  in  Appendix 
B for  comparison  purposes.  Estimated  clockhours  multiplied 
by  team  costs  per  clockhour  established  the  manhour  costs 
for  each  task.  Material  costs  were  determined  by  pricing 
the  materials  consumed  in  each  task.  The  field  level 
tasks  necessary  for  this  analysis  were:  1.  module 

removal  and  replacement/reinstallation,  ?.  module  packing 
for  shipment,  J.  engine  removal  and  replacement/reinstal- 
lation and  4.  base  transportation  to  engine  shop.  Each 
task  will  be  discussed  individually  below. 

Module  removal  and  replacement /reinstallat ion.  Since 
this  thesis  addressed  opportunistic  replacement  at  the 
module  level,  only  removal  and  replioement/reinstallation 
at  the  aggregate  module  level  we re  considered.  Field  level 
module  repair  costs  were  not  considered.  Expendable* 
material  required  for  module  removal  and  replacement  was 
included  to  provide  a total  field  level  module  removal  and 
re  p 1 a c em en t /re  ins  t a 1 1 at  i on  e os  t . 

Module  packing  for  shipment  and  unpacking  upon  receipt,. 
Only  manpower  cost  was  involved  in  packaging  the  module 
for  shipment  once  removed  from  the  engine;  material  is  not 
consumed  during  this  task. 

Engine  removal  and  replacement/re  installation.  The* 
engine  change  operation  essentially  consists  of  three 
operations:  . removal,  0.  reinstallat  ion  and  i,.  trim. 
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Team  composition  cost  multiplied  by  clockhours  required 
determined  manhour  costs.  Interestingly,  an  engine  change 
has  been  made  in  as  little  as  20  minutes,  showing  the 
accomplishments  possible  when  a highly  trained  team  is 
used  and  all  tools  and  equipment  are  prepositioned  (7;9). 

No  consumable  material  was  found  to  be  required  for  engine 
change. 

Base  transportation.  The  time  and  team  size  required 
to  move  the  engine  from  the  aircraft  to  the  engine  shop 
for  teardown  and  repair  constitute  the  basis  for  calcula- 
ting base  transpc  ation  cost.  For  the  present  situation, 
this  cost  is  almost  neglible  but  is  included  for  complete- 
ness. If  a Queen  Bee  operation,  currently  under  discussion, 
is  adopted  for  the  F-100  engine,  base  transportation  costs 
as  defined  would  become  an  inter-base  transportation  cost 
for  the  complete  engine  and  could  be  substantial. 

Module  Failure  Data 

Two  sources  wore  available  for  module  failure 
data — operational  engine  data  obtained  through  the  stan- 
dard engine  status  reporting  system  and  KDf&E  engine  data 
obtained  through  the  F-15  Joint  Test  Force  at  Edwards  Air 
Force  Base,  California.  A complete  computer  program  pack- 
age written  in  FGKTKAN  is  provided  in  Appendix  C to 
manipulate  standard  engine  status  reporting  system  'Affl 
400- ^ ; data.  Acknowledgement  must  be  made  to  the  thesis 
team  of  Pansza  and  Woods  M4j  for  development  of  many  of 
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the  programs  included  herein.  Operational  flying  has  pro- 
duced too  few  engine  failures  to  date  for  any  significant 
results  using  this  reporting  system. 

RDT&E  engine  data  provided  the  only  other  source 
of  engine  removals  for  the  F-100  engine.  A tailor  made 
reporting  system  is  in  use  at  Edwards  Air  Force  Base, 
California  which  provides  data  in  a format  incompatible 
with  AFM  400-1  formating.  Data  was  obtained  in  paper  out- 
put form  from  the  F-15  Joint  Test  Force  YF-100  Engine 
Module  Report.  Data  was  keypunched  onto  punch  cards  for 
input  into  the  CREATE  system,  an  A FLO  Honeywell  6 5b  Dual 
Processor  computer  at  Wright-l'atterson  Air  Force  Base.  A 
final  check  on  data  conversion  from  report  form  to  punch 
card  form  was  accomplished  by  verifying  all  cards.  This 
operation  involves  essentially  retyping  all  data  using  the 
punched  card  just  produced  and  the  original  worksheets. 
Computer  programs  were  developed  to  structure  the  large 
volume  of  data  into  a readable  format  and  screen  it  for 
module  removals  due  only  to  failure  or  expiration  of 
established  operating  hours.  Since  these  programs  were 
developed  for  only  one  time  use, only  the  logic  of  their 
operation  is  presented  in  Appendix  C.  As  future  failure 
data  on  the  F-100  engine  will  be  recorded  in  the  standard 
engine  status  reporting  system;  programs  which  will  screen 
this  data  are  presented  in  some  detail.  A listing  of  KDT&E 
removal  times  upon  which  the  failure  distribution  parameters 
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were  computed  is  presented  in  Appendix  D.  Recalling  the 
discussion  of  potential  hazard  rate  models  from  Chapter 
III,  the  researchers  were  next  faced  with  the  decision 
as  to  which  model  was  most  appropriate  to  model  the 
failure  data  produced.  The  models,  as  recommended  by 
Shooman  are:  1.  exponential  hazard,  2.  piecewise 

linear  hazard  and  J.  Weibull  hazard  (24:194).  Shooman 
further  notes  that  "a  good  way  to  treat  these  (component 
failure)  data  is  to  compute  and  plot  either  the  failure 
density  function  or  the  hazard  rate  as  a function  of 
time  (24:160)."  The  development  of  the  data  hazard  rate 
was  presented  in  Chapter  III  of  this  document.  Shooman 
recommends  the  use  of  "engineering  judgment"  to  select 
the  model  most  appropriate  for  the  data  being  analyzed 
(24:4^7).  Selection  of  the  Weibull  model  was  based  upon 
its  accoptaole  general  fit  when  plotted  against  data 
hazard  computations  and  its  wide  usage  in  reliability 
work  ^24:  *90;56:293  ;.  Furthermore,  the  Wf-ibull  permits 
modeling  both  increasing  and  decreasing  hazard  rates  with 
the  same  mathematical  formula  (24:190).  This  was  of  con- 
siderable value  when  performing  sensitivity  analysis. 

1' lots  of  the  data  failure  rates  are  found  in 
Figures  5* 1 thru  5*4.  The  Weibull  functions  determined 
to  fit  the  data  by  the  program  contained  in  Appendix  E 
are  shown  on  the  same  graphs.  One  notes  the  relative 
"noiseness"  of  engine  failure  data  as  described  by  i'ansza 
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aryl  Voods  previously  (14:45).  This  wide  variation  requires 
additional  caution  when  using  a small  number  of  failure 
points  to  establish  a statistical  distribution  describing 
module  failures. 


CHAPTER  VI 


RESULTS  MD  CONCLUSIONS 


Results 

Initial  data  output.  Figure  H.1  in  Appendix  H is  a sample 
of  the  data  output  provided  by  the  FORTRAN  program 
described  in  Appendix  F.  There  is,  as  anticipated,  a 
point  of  minimum  cost  per  hour.  For  the  baseline  cost 
and  failure  rate  parameters  developed  in  Appendices  B and 
h,  minimum  cost  per  hour  occurred  at  approximately  220 
hours.  The  optimal  cost  per  hour  under  an  opportunistic 
replacement  policy  was  approximately  $.02  less  than  he 
cost  under  a replace  at  failure  policy. 

Figure  6. 'I  is  a plot  of  cost  per  hour  versus  the 
location  of  n for  the  fan  drive  turbine  module  when  the 
replacement  core  module  is  of  age  zero.  Figure  G.2  is  a 
similar  plot  when  the  core  module  is  of  age  240  hours. 

The  effect  of  core  module  age  is  quite  dramatic  wh en  n 
is  located  early  in  the  fan  drive  turbine  overhaul  cycle. 
In  both  the  case  of  a zero  time  core  module  and  a .’40 
hour  old  core  module,  however,  the  curves  are  relatively 
fiat,  in  the  region  near  ?c/)  hours.  Figure  0.2  is  a pLo* 
of  cost  per  hour  versus  the  location  of  n for  the  gearbox, 
given  replacement  core  age  of  zero. 
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Sensitivity  analysis*  The  fan  drive  turbine  nodule  was 
chosen  arbitrarily  for  sensitivity  analysis.  It  is 
anticipated  that  trends  observed  for  the  fan  drive 
turbine  would  also  be  true  for  the  gearbox  module  although 
the  values  of  savings  per  hour,  cycle  length  and  other 
output  variables  would  of  course  be  different. 

Tables  G.l  and  Figure  G.l  in  Appendix  G contain 
a suaaary  of  changes  in  the  location  of  the  optimal  n,  cost 
per  operating  hour,  expected  savings  per  operating  hour 
under  an  opportunistic  policy  and  expected  cycle  length 
as  the  algorithm  input  parameters  were  varied.  The 
results  are  summarized  briefly  in  Figure  6.4. 

For  the  range  of  cost  estimates  and  hazard  para- 
meters explored  when  cost  inputs  or  hazard  parameters 
were  changed  one  at  a time,  expected  savings  under  an 
opportunistic  policy  only  exceeded  $.10  per  operating 
hour  when  the  replacement  core  module  Weibull  m parameter 
was  1.0  or  when  the  replacement  core  module  k parameter 
w*  0.1 0.  With  k = .10,  the  expected  savings  were 
$.1458  per  hour.  With  m = 1.0,  the  expected  savings  were 
$.4317  per  hour.  The  Weibull  m parameter  determines  the 
shape  of  the  distribution  while  the  k parameter  determines 
scale.  In  general,  the  m parameter  primarily  affects 
the  rate  of  change  of  the  hazard  rate  while  the  k para- 
meter for  a given  value  of  m primarily  influences  the 
number  of  failures  per  time  unit.  An  m parameter  value 


Optimal  Savings  Per  K^ur  For  Optimal 

Value  of  n Opportunistic  Policy 
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of  1.0  aod  k parameter  value  of  0.10  were  the  highest 
tested  during  sensitivity  analysis. 

Very  limited  sensitivity  analysis  was  conducted 
varying  two  parameters  at  a time.  We  were  interested  in 
observing  the  effect  on  cost  per  operating  hour  when 
both  variable  costs  were  increased  and  the  core  module 
experienced  either  a relatively  high  m or  k parameter 
value.  With  the  core  module  k parameter  equal  to  0.10 
and  M0D1EEP  (field  level  replacement  cost)  increased 
§100  over  the  base  line  (a  12  percent  increase)  savings 
per  operating  hour  under  an  opportunistic  replacement 
policy  were  approximately  §.17.  When  the  core  module  m 
parameter  was  increased  to  1.0*  and  MODI REP  simultaneously 
increased  by  §100,  saving  per  operating  hour  increased 
to  $.54. 

The  mean  value  of  savings  per  operating  hour, 
for  those  algorithm  input  parameter  combinations  tested, 
was  approximately  §.06  per  operating  hour.  Translating 
this  figure  into  savings  per  year  using  the  fiscal  year 
1981  flying  program  of  175*000  flying  hours  (which  re- 
quires approximately  350,000  operating  hours)  when  F-15 
fleet  acquisition  will  be  complete,  savings  on  the  order 
of  §21,000  per  year  might  be  anticipated  under  an  oppor- 
tunistic maintenance  policy.  We  wish,  however,  to  stress 
that  the  output  data  summarized  in  Appendix  H and  discussed 
in  this  chapter  are  no  better  than  the  hazard  parameter 
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estimates  and  cost  inputs  used.  Hazard  parameter  estimates 
and  coat  estimates  are*  in  our  opinion*  fraught  with 
considerable  uncertainty  at  present.  The  main  interest 
of  t tg  _ * v'vo’v*  of  this  thesis  was  in  developing  a method 
to  find  replacement  policy.  The  data  output 

provided  is  me.  illustrate  only  the  general  range  of 

values  which  l.  Ti . oe  observed.  With  better  hazard  rate 
parameter  estimates  and  better  cost  inputs*  the  relative 
magnitude  of  savings  under  an  opportunistic  replacement 
policy  could  change  significantly.  Illustration  of  this 
possibility  is  found  in  the  behavior  of  cost  per  operating 
hour  when  more  than  one  input  to  the  algorithm  was  changed 
at  one  time.  Further*  there  are  several  inherent  assump- 
tions and  limitations  in  the  algorithm  which  must  be 
understood.  These  assumptions  and  limitations  are  dis- 
cussed in  the  next  section. 

Review  of  Assumptions  and  Limitations 

The  algorithm  developed  in  Chapter  IV  will  deter- 
mine the  optimal  opportunistic  replacement  policy  given 
the  assumptions  and  limitations  which  are  built  into  it. 
The  assumptions  and  limitations  which  were  initially- 
identified  as  pertinent  to  the  research  are  listed  in 
Chapter  II.  Certain  of  these  relate  primarily  to  data 
gathered  in  support  of  the  research.  Those  which  are 
incorporated  into  the  algorithm  are: 
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Module  hazard  rate  functions  were  assumed  inde- 
pendent of  one  another.  This  is  an  aasumption  which 
seldom  holds  true  in  practice.  Certain  types  of  engine 
failures  will  quite  frequently  result  in  damage  to  more 
than  one  component  of  sm  engine.  As  an  example,  foreign 
object  damage  (resulting  from  ingestion  of  a bird  or 
other  solid  object)  may  cause  portions  of  the  inlet  fan 
to  process  through  the  core  module  causing  damage  to  the 
high  speed  compressor.  The  probable  bias  of  the  assump- 
tion of  independent  hazard  rates  would  be  to  increase 
the  number  of  replacement  actions  in  the  interval  (n,  N) 
which  result  from  Joint  failure  of  two  or  more  modules. 

The  convention  used  in  this  thesis  was  to  charge  the  full 
value  of  the  replacement  cost  to  the  module  replaced 
opportunistically  (e.g. , fan  module  or  fan  drive  turbine 
module)  in  the  event  of  this  occurence.  Thus,  relaxing 
this  assumption  (which  would  require  a considerable 
increase  in  the  complexity  of  the  algorithm)  would  probably 
result  in  decreased  savings  per  operating  hour  and  a shift 
of  the  optimal  location  of  n towards  N. 

Transportation  costs  and  depot  costs  were  assumed 
to  be  the  3ame  for  failed  and  unfailed  modules.  As  dis- 
cussed in  Chapter  V,  depot  overhaul  costs  are  not  yet 
well  defined.  A single  cost  for  failed  and  unfailed 
modules  was  used  due  to  lack  of  a more  explicit  cost  infor- 
mation. As  Jorgenson,  McCall  and  Radner  (10:222)  point 
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out*  however,  one  would  normally  expect  the  cost  to  repair 
an  unfailed  module  to  be  lower  since  fewer  components 
should  require  replacement  and  overhaul  actions  should 
be  facilitated. 

A more  subtle  question,  related  to  overhaul  costs, 
was  overlooked  when  developing  the  algorithm.  We  did  not 
include  the  impact  of  changes  in  n on  the  spare  modules 
required  to  fill  the  various  segments  of  the  depot  repair 
and  transportation  pipelines.  Note  that  as  n approaches 
0,  the  expected  cycle  length  would  decrease,  which  would 
cause  required  spare  levels  to  increase.  The  probable 
bias  from  this  omission  would  require  study  to  ascertain 
with  confidence,  for  as  cycle  length  decreased  (and  the 
number  of  modules  depot  overhauled  per  year  increased) 
depot  overhead  costs  would  be  spread  over  a wider  base. 
Given  these  conditions,  use  of  a single  cost  estimate 
for  depot  overhaul  of  a module  would  be  even  more  ques- 
tionable. We  suggest,  however,  without  any  justification 
other  than  intuition,  that  the  cost  due  to  increased 
spares  requirements  would  probably  predominate.  In  this 
case,  the  bias  would  be  to  increase  fixed  costs,  decrease 
savings  per  operating  hour  under  an  opportunistic  policy, 
and  shift  the  optimal  location  of  n towards  N.  Explicit 
consideration  of  the  effects  of  a change  In  n on  spares 
requirements  and  depot  overhead  will,  again,  entail  an 
increase  in  the  complexity  of  the  algorithm. 
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When  developing  the  algorithm*  an  additional 
assumption  which  became  necessary  was  that  all  failures 
in  the  interval  (n*  N)  of  the  module  for  which  opportunistic 
replacement  was  being  studied  result  in  shipment  to  the 
depot.  This  is  probably  not  an  unreasonable  assumption 
for  values  of  n close  to  an  optimally  located  N,  where 
wearout  would  be  significant.  The  assumption  is  essen- 
tially indefensible  as  n approaches  0,  however.  If  the 
assumption  were  valid*  there  would  be  scant  justification 
for  an  intermediate  level  module  repair  capability.  We 
suggest,  however,  based  on  results  of  the  sensitivity 
analysis,  that  the  region  of  interest  where  the  optimal 
value  of  n is  most  likely  to  be  found  is  near  N rather 
than  near  0.  The  probable  bias  resulting  from  the  assump- 
tion that  module  failures  in  the  interval  (n,  N)  are 
depot  reparable  only  would  be  to  increase  expected  cycle 
length  under  an  opportunistic  policy  and  decrease  the 
probability  of  a field  level  replacement  action.  Both  of 
these  effects  would  result  in  a decrease  in  cost  per 
operating  hour  under  an  opportunistic  replacement  policy. 

There  is  an  inherent  limitation  in  solving  for  n 
independent  of  N.  To  illustrate,  during  sensitivity 
analysis  the  maximum  operating  times  for  the  fan  drive 
turbine  and  core  modules  were  experimentally  increased  to 
500  hours.  When  this  change  was  incorporated,  optimal 
cost  per  operating  hour  dropped  to  $56.91  which  is  eight 


79 


percent  leas  than  that  achieved  with  HOT*  a equal  to  230 
hours.  Baaed  on  the  results  o f the  study  by  Jorgenson* 
McCall  and  Radnor  (10:223-235),  it  ia  economically  optimal 
to  solve  for  n and  N simultaneously.  Only  a moderate 
change  to  the  FORTRAN  program  in  Appendix  F — including  an 
additional  do-loop  to  scan  over  the  range  of  interest  for 
maximum  operating  time  and  additional  statements  to  store 
minimum  costs  as  the  do-loops  were  executed— would  be 
required,  but  the  hazard  function  must  be  clearly  defined 

•i. 

over  the  relevant  region. 

At  the  time  interviews  were  conducted  to  gather 
cost  data  (20),  the  subject  of  test  cell  procedures  for 
the  F-100  engine  was  not  well  settled.  The  original  con- 
cept was  to  accomplish  engine  trouble  shooting  on  the 
aircraft  prior  to  in  shop  maintenance.  There  is  at  the 
present  time  a trend  towards  engine  trouble  shooting  on 
the  test  cell  prior  to  an  engine  undergoing  maintenence. 
Test  cell  cost  was  not  included  in  the  calculation  of 
base  level  replacement  cost.  One  would  expect  inclusion 
of  this  cost  to  result  in  increased  savings  per  hour  under 
an  opportunistic  replacement  policy,  since  it  is  a variable 
cost,  and  movement  of  the  optimal  value  of  n away  from  N. 

Conclusions 

With  the  limitations  and  assumptions  outlined 
above,  an  algorithm  has  been  developed  to  determine  the 
optimal  location  of  the  breakpoint  between  the  replace  at 
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failure  region  and  the  opportunistic  replacement  region 
of  the  module  overhaul  to  overhaul  cycle.  This  answers 
research  question  number  one.  Sensitivity  of  the  optimal 
location  of  the  breakpoint  to  changes  in  hazard  rate 
parameters  and  changes  in  ccst  inputs  was  examined  through 
sensitivity  analysis.  The  optimal  location  of  the  break- 
point varied  from  the  maximum  operating  time  to  30  hours 
less  than  the  maximum  operating  time  for  those  values  of 
input  hazard  rate  parameters  end  cost  inputs  studied. 

This  answers  research  question  number  two.  As  the  inputs 
to  the  algorithm  were  varied,  expected  savings  possible 
under  an  opportunistic  policy  varied  from  less  than  8.01 
per  operating  hour  to  slightly  over  8.30  per  operating 
hour.  Only  in  the  case  of  a core  module  Weibull  m para- 
meter equal  to  1.0  did  expected  savings  exceed  8.13  per 
hour.  Thus,  expected  savings  per  hour  under  an  oppor- 
tunistic policy  is  relatively  insensitive  to  changes  in 
cost  and  hazard  rate  parameter  inputs.  This  answers  research 
question  number  three.  The  average  magnitude  of  expected 
savings  under  an  opportunistic  policy  was  8.00  per  opera- 
ting hour  for  those  hazard  rate  parameters  and  cost  inputs 
studied.  8.00  per  operating  hour  represents  an  approxi- 
rate  0.1  percent  savings  over  a replace  at  failure  policy. 

In  terms  of  the  FY-81,  F— 1 5 flying  orogram,  savings  on 
the  order  of  821,000  might  be  expected.  This  answers 
research  question  number  four. 
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Re  c ommenda t i ons 

Although  the  authors  would  like  to  believe  that 
this  thesis  reflects  a moderate  amount  of  research  effort) 
we  recognise  that  we  have  barely  scratched  the  proverbial 
surface  in  terms  of  the  amount  of  work  yet  to  be  done  in 
the  area  of  aircraft  engine  module  opportunistic  replace- 
ment policy.  In  particular*  we  would  suggest  that  effort 
be  directed  toward:  1.  improving  hazard  rate  parameter 

estimates  and  cost  estimates  as  the  P-15  aircraft  and  P-100 
engine  accrue  more  operational  experience,  2.  incorporating 
necessary  changes  into  the  algorithm  so  that  some  of  the 
restrictive  assumptions  employed  in  this  thesis  can  be 
relaxed,  3.  exploring  the  effect  of  changes  in  cycle 
length  on  spares  requirements  and  depot  overhead  charges 
and  4.  performing  additional  sensitivity  analysis  by 
varying  more  than  one  input  variable  at  a time  and  expan- 
ding the  range  within  which  variables  are  tested. 
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DEFINITIONS 


Augmented  Turbofan  Etagine  - A basic  turbofan  jet  engine 
with,  an  augment  or  attached  to  discharge  end  of 
turbine  section.  An  augsentor  nixes  the  hot  turbine 
discharge  gases  and  the  relatively  cool  fan  bypass 
air.  The  mixture  thus  obtained  is  burned  in  the 
after-burner  segsents. 

Catastrophic  Failure  - A failure  characterized  by  sudden * 
unexpected  daasge  or  loss. 

Concurrent  Maintenance  - Accomplishment  of  two  or  more 
independent  maintenance  actions  at  the  same  time. 

Cycle  Time  - The  length  of  time  from  installation  of  a 
module  until  its  removal  for  maximum  operating  time 
(MOT)*  failure*  or  opportunistic  replacement. 

Deterministic  - The  result  of  a given  action  is  known 
with  complete  certainty. 

Downtime  - Any  time  period  in  which  an  aircraft  or  com- 
ponent is  not  available  for  use.  This  is  normally 
classed  as  Not  Operationally  Ready-Maintenance 
(NORM)  or  Not  Operationnaly  Ready-Supply  (NORS). 

Echelons  - Levels  of  the  maintenance  organizational 
hierarchy. 

End  Item  - An  item  selected  for  specific  configuration 
and  accounting  control,  (e.g. , aircraft  engines 
and  airframes;. 

Fixed  Costs  - Those  costs  known  to  occur  with  certainty 
during  a single  cycle.  For  this  research*  fixed 
cost  included  engine  removal  and  installation  costs* 
module  packing  and  unpacking  costs*  and  transporta- 
tion costs  for  a module  between  bas'  and  overhaul 
facility. 

Lost  Service  Life  - The  sacrificed*  otherwise  available* 
service  life  of  a component  due  to  replacement  before 
failure. 


85 


84- 


Module  Construction  - The  concept  where  functionally  and 
physically  associated  parts  are  removable  as  units. 

Opportunistic  Replacement  - The  replacement  of  an  end  item 
specifically  during  a time  when  the  item  must  be 
removed  to  perform  other  maintenance  (e.g. , replace- 
ment of  the  removed  inlet/fan  module  given  mandatory 
core  replacement). 

^ueen  Bee  Operation  - The  maintenance  concept  of  assigning 
intermediate  level  maintenance  capability  for  specific 
r-ype,  model  and  series  of  aircraft  engine  to  a 
centralized  location  which,  in  turn,  provides  support 
to  other  bases.  Under  this  concept,  supported  bases 
would  basically  have  only  remove  and  replace  or  minor 
repair  capability. 

stochastic  Failure  - The  time  at  which  failure  occurs  is 
not  known  with  certainty,  i.e.,  the  time  at  failure 
is  governed  by  a probabilistic  mechanism.  Only  the 
expected  time  of  failure  can  be  determined. 

Test  Cell  Kun  - Operation  of  a jet  engine  on  a specially 
designed  fixture  (”0011;  after  removal  from  an  air- 
craft fcu*-  prior  to  teardown,  specifically  accomplished 
to  isolate  a defective  component  or  components.  Test 
cell  runs  also  are  made  after  build-up  or  repair  to 
verify  maintenance  actions. 

Trouble-shoot  - The  maintenance  actions  necessary  to 
isolate  a defective  component  or  components.  This 
may  be  accomplished  on  the  aircraft  or  after  removal. 

Variable  ,'osts  - Costs  which  vary  in  amount  or  may  or  may 
not  occur  during  a single  cycle.  As  an  example,  in 
the  event  of  opportunistic  replacement,  certain  costs 
are  nor.  incurred  which  would  be  incurred  in  the  event 
of  a replacement  at  failure. 

Wearout  - The  characteristic  where  hazard  rate  increases 
with  age. 
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Table  B.c 

Standard  Pour  Man  Team  Coats 


Wage  iiate 

Personnel  Cost 
Ciockhour 

84.75 

8 4.7^ 

85.9 

8 3.92 

85.29 

8 6.58 

8 '5.25 

i er  hour  worked.  Source:  AHYi  177-"  01 


Tab!*-  H.5 

standard  it.ree  Han  Team  Costs 


Hank 

Wage  Hate 

Cgt 

A :'J 

85.92 

85.29 

Personnel  Cost  Per 
Clockuour 


* 5.92 
8 6.58 


>tal  Team  Cost /Hour . 8/'0.50 


Per  Hour  worked.  Source:  AFH  1 77-^01 
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Table  B.4 

standard  Two  Kail  Team  Costa 


Number  Required 

Rank 

"1 

Wage  Rate 

Personnel  Cost  Per 
Clockhour 

Jt 

\ 

Sgt 

*5.92 

*5.92 

1 

A1C 

*5.29 

*5.29 

Total  Team  Cost/Hour *7.2' 


* 

'Per  hour  worked.  Source:  AFP!  177-*01 


Replaceaent/Reinstallation  Material  Coats 


Jnit  Price  I Total  Price 
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1 Estimates  were  15-30  minutes  so  the  average  was  used. 
Total  RScR  Engine  Costs  - 1161.84-  (which  is  the  sum  of  total 
costs  from  Table  8 and  9). 
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APPENDIX  C 


PROGRAMS  USED  TO  SCREEN 
MODULE  FAILURE  DA  *1*1 


Apparoxx  c 


ISOGRAMS  USED  TO  SCHEES 
MODULE  FAILURE  DATA 

This  appendix  contains  a copy  of  each  different 
type  of  computer  program  used  in  the  module  failure  data 
screening  portion  of  the  research  effort.  The  programs 
included  here  are  representative  of  the  programs  used; 
i.e.*  some  programs  used  were  duplicates  of  the  programs 
here  with  the  exception  of  the  data  files  processed. 

The  manner  of  presentation  used  for  this  appendix 
will  be  to  list  the  variables  and/or  files  used  in  each 
program*  followed  by  an  actual  listing  of  the  programs. 
Listing  the  programs  was  made  possible  by  use  of  a 
computer  program*  NICELIST*  developed  by  Major  Jim  Abbott « 
Computer  Support  Section*  School  of  Systems  and  Logistics* 
whom  the  authors  sincerely  thank. 

Many  of  the  computer  programs  listed  in  this 
appendix  and  used  in  the  research  effort  were  developed 
by  Fanasa  and  Woods  (14:95-^21 ) and  are  used  with  their 
permission. 
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VOTE:  This  CARBIg  Pro^jy^canrerta  data  tram  BCD  tape  to 

UEATT  — The  sane  tic  tape  containing  the  trans- 
action history  of  all  operational  F-100  engines  through 

3lft  .Tima  ~ 

50  «iuns  T5r/?« 

HBftAFE  - The  permanent  disk  file  containing  the 
data  generated  by  the  progrea  TAPEPILE.  This  data  was 
obtained  ffoa  the  nagaetic  tape  VTATT. 


Program  TAPEFILE  Listi 


81. 


om  Mtir*  data  vMoid  into 


A1  - Nodal*  Serial 
A2  - Station  (Boo*) 

A3  - Bat*  of  Tin**rt1nn 
At  - Nodal*  Transaction  *ad 


Condition, 
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A5  - Nodal*  Bonotal  H***on  aad  Nodal*  Hour*  Sine* 
Overhaul. 

A6  - Bogin*  Designation  and  Ehgine  Serial  Buaber. 

H - A counter  used  to  indioat*  the  nuaber  of  record* 
processed. 
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A - Variable  in  res ding  the  fint  two  da ta 

elements  of  the  data  fila  Wlfflrtg.  Tteee  data  elements 
have  no  site** icance  in  this  research  effort  but  are 
retained  to  simplify  formatting  compatibility  with  otter 
proftaH* 

TTPWOD  - Variable  used  to  specify  to  which  of  the 
five  modules  this  data  element  refers.  VA  refers  to  a fan 
module  entry.  JB  refers  to  a core  module  entry.  VC  refers 
to  a fan  drive  turbine  module.  VD  refers  to  an 
augmaator/cteauat  module.  KE  refers  to  a gearbox  module. 

69  refers  to  an  entire  engine  entry. 

B - Variable  used  in  reading  the  last  45  data  ele- 
ments in  the  data  file  HEWTAPE2.  These  data  elements 
contain  information  such  as  data  of  report,  reason  for 
report,  reason  for  removal  (if  appropriate),  engine 
operating  hours  as  of  this  date,  engine  serial  number. 

I - A counter  used  to  indicate  the  total  number  of 
data  records  processed. 

J - A counter  used  to  indicate  the  number  of  inlet 
fan  module  records  processed. 

K - A counter  used  to  indicate  the  number  of  core 
module  records  processed. 

L - A counter  used  to  indicate  the  number  of  fan 
drive/turbine  module  records  processed. 

M - A counter  used  to  indicate  the  number  of 
augment or/exhaust  modules  records  processed. 

N - A counter  used  to  indicate  the  number  of  gear- 
box module  records  processed. 

NEETG  - A counter  used  to  indicate  the  number  of  whole 
engine  records  processed. 

II  - Variable  used  to  sun  the  number  of  records 
written  to  module  files  ted  the  engine  file. 
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Di.  Pronraa  SSQLST  Variables  asd/"  Tlltt  7ft1 

A - The  file  line  number  associated  with  one  record. 

E - The  Module  Serial  Stihber.  In  this  format  the 
naaufaetor’s  code  and  all  leading  zeros  have  been  strip- 
ped off. 

V - Defined  all  the  data  in  one  record  except  that 
defined  by  the  variables  A and  E above. 

SEE  - Module  serial  nuaber.  A mechanism  used  to 
d&t ermine  if  the  transaction  just  read  by  the  program 
SHOLST  was  related  to  a Module  Serial  Suaber  that  was 
different  from  the  one  previously  read. 

Z - A counter  used  to  indicate  the  nuaber  of  records 
processed. 

OCAMAl  /DATA1  /XH?AV  - ASC  II  permanent  file  structure 
containing  the  transaction  history  of  the  inlet/fan 
modules.  OCAMAl  is  a catalogue  of  all  data  from  the 
D024  Ehgine  Status  Reporting  System.  DATA1  is  a sub- 
catalogue containing  all  data  prior  to  screening  for 
errors.  UFA*  is  the  file  asms. 


El.  Program  ZEBOPT  Variables  and/or  Files  Used. 

A - The  file  line  number  associated  with  one  data 
record. 

B - Module  Serial  Humber. 

C - Engl ne  designator  and  Module  Removal  Reason. 

D - Module  Operating  Time. 

E - Defined  all  the  data  in  one  record  except  that 
defined  by  the  variables  A,  B«  C and  D above. 

SEN  - Module  Serial  Number.  A mechanism  used  to 
determine  if  the  transaction  just  read  by  the  program 
ZEROPT  was  related  to  a Module  Serial  Number  that  was 
different  from  the  one  previously  read. 

N - A counter  used  to  indicate  the  number  of  records 
processed. 

OCAMA1/DATA1  /INPAN  - ASC  II  file:  structure  as 

defined  in  Program  SNOLST. 


E 2.  Program  ZEROPT  Listing 


-J 

< V2 

t-4  t— < 

* 

2 :r: 

Ul  H 

0 

«1 

71  1-4 

• 

**► 

— 

H _ 

• 

* — * 

to  c 

s 

• 

oc  C 

00 

z 

t-4  >1 

**• 

to 

UL  Ct 

< 

< 

LU 

lu 

Ul  Q. 

2 

rc  'O 

t-4 

Q 

>4  iu  * 

V 

Ul 

<— ■ 

to 

(/)*->  LU 

< 

HH  (O 

H 

71  • 

< LU 

4"4 

C 

*-•  u;  i" 

C * 

r* 

LJ 

-J  o :u  cm 

V,  C 

»v 

•-* 

-*o  * 

c 

u 

« H2  Ul 

~<  rj 

0 

C 

Ul 

H < Ul  * 

r-  ?♦  > 

V 

IT 

u; 

C\<  uL 

< 

c ai  Cl  + 

<o  • 

0 

W» 

LU 

rxauic 

LU  O 

H 

tu  ci  to  * 

V S « 

0 0f. 

•Nl  O 

<IC  o. 

«* 

c • 

* 

vSX« 

*n  <S 

* 

• o 

3 

ill 

a ui  o 

«f>  | 

**■% 

o < 

r^. 

? mx  * 

<—  n 

CO  O 

• • 

S 

£ 

< K X.  O 

► S.  2 

_ H 

•n  >0 

-** 

r» 

ne  2 * 

O:  lu 

«W* 

— 4ft 

•4| 

0 

**♦ 

0 w-  «•—  £— * 

< » * 

H w 

C • 

»> 

c u * 

•—  V 

a 

• X 

s£ 

lU 

4<  JO 

>c  — 

^ **** 

8 

*-i  << 

tn  <*•  re 
to  ui  a 

t—  to  uj 

m h 

isra: 

H-5  H *<  ■: 

5£ 


P5 

<< 


— u.l 

— (0 

• • 

o c 

...  „ Ml  LU 

c<  * * • 

u.  !u—  mn 
■—I  a:  ♦ w *-> 
*■?  _i  *r  u.  ul 

■■2  C <*  tl  *■*  w* 

— O 


"IkH 


— a 


CM  «-h 

«c 

au  f~  t: 
« 2 c: 
2 •—  c 

? * I £T  U. 

•, o o. 

^ ^ 
ro  n cv 


Ul  m 
— . •*>  U 

— n—  £ 

W 1-4 

ottcrc 

H *-•  H ft. 


- ..  ..  - *5  s 


>•*»  «<■ 

CM  f'  V 


* 


© © 5J  © 

ITS  Of- T 


CL 

c a 

H* 

to  LU 


o *r- 

— IV 


112 


Pi . Program  BRCHBCK  Variables  and/or  Files  Used. 

A-  Line  Number , Module  Serial  lumber  and  Module 
Designator. 

B - Module  Removal  Reason. 

C - Defined  all  the  data  in  one  record  except  that 
defined  by  the  variables  A and  B above. 

RM  - A vector  used  to  store  each  different  module 
removal  reason  during  the  execution  of  that  part  of  the 
program  BRCHECK  which  identified  these  values.  Also 
used  to  provide  a listing  of  the  different  values. 

N - A counter  used  to  indicate  the  number  of  records 
processed. 

J - A counter  used  to  index  the  vector  RN. 

RPAIRC  - ASC  II  permanent  file  containing  the  trans- 
action history  of  those  sample  modules. 
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oi.  ’Trana  7™nr”  trf /o-  tUa 

B - The  file  line  number  associated  with  one  record. 
C - Type  module  code. 

0 - Zeroes. 

£ - Module  Serial  Huuber 
F - Station  (Base)  Haae. 

6 - Bate  of  Transaction. 

H - Module  Transaction  and  Module  Condition. 

J - Module  Beaoval  Season. 

K - Zero. 

M - Module  Hours  Since  Overhaul. 

N - Blank. 

P - Engine  Type. 

S and  T - Engine  Serial  Number. 

Z - A counter  used  to  indicate  the  number  of  records 
processed. 


ZEROCKN  Listi 


A - Line  nillur  associated  with  one  file  record. 

B - Hodule  Serial  lumber. 

BB  - Hodule  Desipiator. 

C - Hodule  Removal  Beeson. 

0 - Hodule  Operating  Tine. 

E • Defined  so  ell  the  date  on  one  record  except 
those  defined  by  the  variables  A*  B«  BB*  C*  and  D above. 

SV  - Hodule  Serial  Humber.  ▲ nechanisa  used  to 
determine  when  all  of  the  transactions  pertaining  to  a 
given  Hodule  Serial  Humber  have  been  processed. 

OFT  - Hodule  Operating  Tine.  A nechanisa  used  to 
detect  an  erroneous  decrease  in  operating  tine. 

V • A counter  used  to  indicate  the  total  number  of 
records  processed. 

WYATT  - ASC  XI  permanent  disk  file  that  contained  the 
screened  module  transaction  histories. 
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II. 


A - Lin* 
Designator. 


lift  and/or  Files  Used. 

• Module  Serial  Mudber  and  the  Module 


B - Module  Henoval  Beeson. 

C - Module  Operating  l.jie. 

D - Defined  as  all  the  data  eleaeata  in  one  record 
except  those  defined  in  the  variables  A,  B end  C above. 

B • A counter  used  to  indicate  the  t&tel  of 

records  processed. 

WXATT  - Permanent  disk  file  that  contained  the  trans- 
action history  records  of  all  nodules. 


t 

a 


i 


i 


i 

i 


1 


7 

i 

I 


s 


120 


Jl.  Pto ar—  DUPCA 

A - Line  lumber,  Nodule  Serial  lumber  end  Nodule 
Designator. 

B - Nodule  Removal  Reason. 

G - Nodule  Ope refine  Tine. 

U - Defined  as  all  the  data  in  one  record  except 
that  defined  by  the  variables  A*  B and  C above. 

V - A counter  used  to  indicate  the  total  number  of 
records  processed. 

H - U*>d  as  an  indication  to  the  program*  that  the 
first  transaction  with  a removal  reason  had  been  processed. 

TIME  - Nodule  Operating  Time.  A mechanism  used  to  detect 
duplication  of  module  operation  time  points  between  suc- 
cessive module  removal  transaction. 

WYATT  - Permanent  disk  file  that  contained  the  trans- 
action history  records  of  all  modules. 
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K1.  Program  REMOVALS  Variables  and/or  Files  Used. 

FILENAME  - Variable  name  used  for  the  name  of  the 
data  file  to  be  read  froa. 

RITENAME  - Variable  naae  used  for  the  naae  of  the 
data  file  to  be  written  to. 

A - Defined  all  the  data  in  one  record  for  the  first 
23  spaces. 

REMOVAL  - Variable  used  to  contain  Removal  Code 
entries  from  data  file. 

B - Defined  all  the  data  in  one  record  for  the  last 
23  spaces. 

I - Counter  used  to  count  nuaber  of  records  read 
from  FILENAME. 

J - Counter  used  to  count  nuaber  of  records  written 
to  RITEVXLE. 


123 


5 

4 J 
HIS 

uizPa. 

Bn 

JK  LUX 
LUCBHUI 

^tSS5 


XLUUJX  * 
tx  XHuj 
N p<  s 

in  x ec< 

JOB  ^ 


LU 

U 

x 

< 

(2 

>•«*• 

dd 

g£ 

cs 

I IS 

i 

S'* 

ts 

t-4 

Sfe 

5“ 

a.  lu 


2 

> 

Is 

x^* 

5 UJ  • 

atn 

g"ffl 

coB 

tux  < 
h 

< x cs 
UJmUJ 

a a 

O V. 

• a 
>SB 

uihSc  s 

SnnUJ  •» 

a * » * 
< 3 o • x 
in  < 

IUUICmU, 
>®3-  2S 

< < H*- 

xc  P V 
huSlA 
P OBJ 
CV)2n< 
>•*-..  > 


>LU  Jh 
O UJNN 

IhJlLS 

CD 

• 

Bfes 

S •<■»  s** 

P 

x<  < o 
ahiffl 
< < 5 a 

xp  m • 

* 

a 

x.  r?£  a 

. _ ® LU 

V JUi * > 

§ 

225*3  $ 

See  <zn  • 

a in  ui  « n 

2c~i2 

XX  _ - _ 

p 

as  ace  z 

UJUitUM 

sm 


mo.  ex 
os  c 


H 


tn  r\ 

ESS* 


tn  < 

UJ.T 

Qp 

pa 

5 uj 


<n  »j 
• Z H < 


a.  > 
Ifl«sc 

» 


x.§3£ 


- < - 

a aa  _ 

uub« 

K^i 


ia 
fc££ 

•n  ck  a 

a sun. 

**■  p © 

*9  *5  — 

» — — 


ID 

Uj 

>* 

a. 

x 

3 

in 


x a c — 
j jao< 
•-•  < x 
x > in  x < 


tn 

x 

>• 


:mfS  uj\ 

HaS^ 


:u 

<W2 

*•§  in 

*>*  «e  <■* 

a iu  x 

a.  * >m 

® 

<v  n «♦ 


iu*S!C* 



uj  in  h:  ii  o < 
a.  •->  o < — x 
>•  x an  x a.  :u 
HH<XSH 

> • » • » • rr. 

fefefefet^Q' 


IU 

p 


« 

X. 

H 


s> 

Urn  X 

M 

• tn 

X 
ui  X 

X 

•.u 

x a 

2o 

c5° 

J IU 
— C 

*fcuj 

X X 
p a 

§< 
X X 

< UJ 

t 1 

ox 

w 

a;  x 
a.  c: 
HO 

rr  ~ 


3e 

x 


< 

c 

\ 


3 

o 

\ 

X x 
n z 
a < 


btn  — r*. 

ft  Iff 


a a re  ns  nc  mu 
a.  n.  n.  o.  a o.  rr 

a s,  c ts  3 o o. 

incMBo*!.  — 

— — — — — fvi  f\ 


N*«  llj< 

H s * * 
Z » • » 


O 


f»  rr.  rr  nc 
l%  0»  P CL 


ji  sv  ^ 

f VPItr/l 
PV!  PVJ  fV  IM 


a.  x' 
\j: 

• X *~) 

H • •* 

ft-ft  • * 

ti’ 

rcuj-* 
n rv  *-» 

S*  ~ ^ 

« fx*  T) 
<V  PV  IV 


V 

w 

H 

S.v.: 

r 
H **■ 

•**  iL. 

a « 
a tc 


«»  O* 

o i? 

nj  m 


3U*  CALL  A7TACF(!3tFILENAMEt3.0.t 
32<>  IFC  ANS*NE«MYES,,)G()  T > 24 
33C  CALL  ATTACH(l,HITENAVft:,3#a,#> 
34tf  CALL  FMcDI AC  I ,5) 


♦ 

Ui 

> «— 
S f\l  o 

• — c 
< 

—»  • 

# a 
C < u 

0 2*3 

■ U4  X • 

5> 


?/■)  *- 
■ 2 UJ  C 

;<Hh 

lifts 

l“T(5 


a 

a 

t/>  o- 

*■<  c: 

UJ 

a cr 

• “a  •— 

uj  »r> 

~~  a O 
< LU 
ZWK 
UJ  M 

u (/) 

• <<  *M 

a »(/> 

ziS 

</>  c o 

Q UJ  V O 

tt  J 

C Z7U  o 
U C < H 
tilOU 

«**  ML’  •* **  #-■» 

- -r,  v 

•-4  0 

LU  XH  — 

w ^ ^ 5 

Q t/7 

iXJU  UJ 
[LI  UJ  < > 

llj  LL  S 

r:  ?:  m • 

1 D3tr  c 

7*  * 2 ."5  < UI 

1 • * s s a a a * 

'►—•••••<  rn 

JHH-HH  *2 

I**-*— 


fva.fl.flLtt.CLa?— »v?uj 


ssi&osca&iSjiSi 


3 3 Ui  U CSJ  ca 


* 4 


I » 


125 

Programs  to  handle  Edwards  Air  Force  Bfus*  P,')TfeE 
Engine  Data  were  developed  to  manipulate  the  data  provided 
into  the  proper  format  for  usage  beginning  with  Program 
SPLIT-OC . The  programs  listed  in  this  appendix  were  then 
used  to  screen  all.  files  of  obvious  errors.  Since  the 
initial  manipulation  programs  were  designed  for  one  time 
usage  only,  the  authors  have  not  included  them  in  this 
document. 


APPfUDIX  D 


MODULE  TIMES  AT  FAILURE 


Number  of 
MOT  Rmvls 
Recorded 


1.6 

2.9 

6.1 

9.8 

18.2 

23.9 

27.6 

29.5 

32.6 

32.7 
49.0 

101.8 

112.1 

122.0 

133-6 

159.2 


14.6 

21.7 

64.0 

60.8 

86.7 

127.1 

180.5 

164.7 

237.5 

174.4 

207.4 

220.4 

220.7 

2 

8 

18 

fl 

Total  No. 
of  DATA 
Po  ints 


14 


r 


m Parameter 
estimate 

Std  Error 

'l 

of  Estimate 

ggSggSBS S8SS9ESB3B 

Parameter 

'l 

Estimate 

Hfcd  Error 
| of  Estimate 


-0.2407 

-0.0132 

0.6601 

0.71019 

0.0829 

1.1106 

0.1250 

0.0800 

| ggg| 

0.00205 

0.00014 

0.000163 

0.0070 

0.00092 

0.00005 

. 

0.000043 

As  determined  using  LIKELY2,  Appendix  E. 
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APPENDIX  £ 


DEVELOPMENT  OP  A COMPUTES  PROGRAM  TO  PERFORM  MAXIMUM 
LIKELIHOOD  ESTIMATES  (MLE)  OP  THE  PARAMETERS 
OF  A TRUNCATED  WEIBULL  DISTRIBUTION 


apphtdix  e 


DEVELOFtlflfT  Of  A COMPUTES  MOGBAM  TO  PERFORM  MAXIMUM 
LIKELIHOOD  ESTIMATES  (MLE)  OF  THE  PARAMETERS 
OF  A TRUNCATED  WEIBULL  DISTRIBUTION 

In  working  with  data  from  component  life  cycle 
testing*  it  is  usually  desirable  to  generalize  from  sample 
data  to  the  entire  population  of  similar  components. 

Shooman  (24:195)  points  out  the  essentiality  of  fitting 
failure  data  to  a statistical  distribution  for  this  purpose. 
Shooman  further  urges  use  of  maximum  likelihood  estimators 
for  parameter  estimation*  once  a specific  distribution  is 
selected,  since  they  offer  "the  most  flexible  and  powerful 
of  modern  estimation  techniques  (24:472)."  As  discussed 
in  Chapter  V*  the  authors  chose  the  Weibull  distribution 
for  this  research  effort.  The  specific  form  of  the  Weibull 
distribution  used  is: 


Probability 
Density  : f(t) 

Function 


Kttte  “ Kttt+1/(m^) 


f t,' 

vh.  ■ j 


m >-1 

K i 0 


where  m and  K are  the  parameters  of  the  distribution  and 
t is  the  variable  of  interest— in  this  situation— time. 
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Once  failure  data  wee  screened  and  data  points 
established  (Bee  Appendix  C for  procedure  used)*  the  task 
became  to  determine  the  appropriate  parameters  for  the 
Veibull  distribution  to  nodel  each  nodule  hazard  rate. 

One  notes  the  existence  of  an  established  maximum  opera- 
ting tiae  (HOT)  for  each  module*  except  the  augment or/exhaust 
nodule.  This  HOT  essentially  truncates  the  distribution 
and  this  effect  aust  be  addressed  in  order  to  have  valid 
parameter  estiaates. 

Maximum  likelihood  estimates  (MLE)  of  the  ehape 
and  scale  factors  of  a truncated*  Weibull  distribution 
have  been  developed  by  Shoonan  (2<f:4 77)*  The  expressions 
for  the  shape  factor*  n*  and  scale  factor*  X*  where  r is 
the  number  of  failures  in  tine  T*  n is  the  nunber  of 
components  in  the  original  population  and  t^*  t2»  • • •* 

t are  the  failure  times  of  the  r components  are: 

K - ■■■  -Stl+U (E.2) 

£ t^1  ♦ (n-r)Ttt+1 

i-1 


Lln  ti 

i-1 


r 

E 

i-1 


*i 


m+1 


In  t,  + (n-r)!**1  InT 


£ tf*1  ♦ (n-r)!**1 
i-1 


1 

r 


1 

mTT 


(E.5) 
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Equation  (E.3)  cannot  be  solved  explicitly  Tor  a. 
Shooman  (24:477)  recommends  graphical  solution , that  is 
plotting  the  left  side  of  Equation  (E.3)  against  the  right 
side.  Although  this  is  certainly  an  alternative,  the  cur- 
rent authors  decided  to  take  advantage  of  the  computational 
power  of  the  computer  and  develop  a PORTRAIT  program  which 
would  search  for  a value  of  m which  would  satisfy  Equation 
(E.3).  It  can  be  noted  that  the  left  side  of  Equation  (E.3), 
not  containing  m,  is  a constant.  This  simplified  the 
programming  in  that  only  the  right  side  of  the  equation 
needed  to  be  repetitively  recalculated.  A graphical  chart 
depicting  the  computer  search  operation  is  included  in 
Figures  E.1  and  E.2.  A copy  of  the  program,  LIKELY 2 is 
included  later  in  this  appendix.  In  order  to  verify  that 
the  program  functioned  properly,  a number  of  data  files 
with  known  Weibull  shape  and  scale  parameters  were  created 
using  the  inverse  transform  method.  One  file  of  1000  data 
elements  was  created  with  a positive  shape  parameter  and 
one  of  the  same  a ize  with  a negative  shape  parameter. 

With  n set  equal  to  r and  the  truncation  point  set  Just 
above  the  value  of  the  largest  data  element,  LIKELY?  was 
used  to  estimate  the  parameters  of  the  underlying  distri- 
bution. Finally,  the  positive  shape  parameter  file  was 
truncated  by  simply  splitting  the  file  into  two  smaller 
files,  one  containing  values  greater  than  the  truncating 
point  and  one  smaller  than  the  truncatinp  point.  ILKELY2 


ALT  (see  Page  1J6) 


Figure  E.1 

Graphical  Depiction  of  LIKELY2  for  Positive  a Search 
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than  used  to  estimate  pmaattn  from  tbe  truncated 
data  file.  Basalts  returned  tQT  LTHT»T2  were  found  to  be 
statist  leal  ljr  accurate  at  tbs  95  pare  sat  confidence  level 
when  compared  to  tbe  paraaatars  usad  originally  to  writs 
tha  data  points. 

formulae  ara  also  available  (34:482-48?)  to  find 
tba  amount  of  variance  involved  in  each  estimate  of  tba 
paraaatars  a and  K.  Tba  expressions  for  parameter  esti- 
mation variance  for  a truncated  tfeibull  ha  card  model  where 
a is  tba  estimate  of  tba  shape  parameter*  K is  the  estimate 
of  tbe  scale  parameter*  n is  tbe  number  oi  components  in 
the  original  population*  r is  tbe  number  of  observed 
failures  in  tbe  test  period  T*  and  tA*  tg*  . . . » tr  are 

tha  failure  times  for  tbe  r components  are: 

Var  K n ^/r  (E.4) 


Var  a * 


1/K 


] 


(E.5) 


0+1 


la  tt  + (n-r)Ttt‘*‘1  InT 


] 


H«i 


i®*1  (lnt^2  ♦ 


t m ^ jm+i  (3^)2 


i-i 
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Using  Chebyshev' s inequality  (24:483)  which  states 

that: 

P(|#-s|*N#)*p  (£.6) 

where: 

6 random  variable  with  an  arbitrezy 
distribution 

a expected  value  of  6 or  E(  9 ) 

N number  of  standard  deviations 
standard  deviation  of  9 or  ^ 
one  may  determine  appropriate  confidence  bands  for  the 
parameter  estimates. 

The  computer  program  developed  in  FORTRAN  to  find 
MLE’s  and  their  variances  is  shown  in  Figure  £.4.  Basically 
the  program  consists  of  a Main  Program  and  seven  function 
subprograms.  The  main  program  initializes  all  variables, 
reads  data  elements  in,  determines  when  m has  been  esti- 
mated to  an  accuracy  of  .00001  and  provides  output. 
Inspection  of  the  right  hand  side  of  Equation  (£.33 
revealed  its  monotonic  behavior.  Because  of  this  behavior, 
the  interval  bisection  offered  considerable  computational 
efficiency  and  was  incorporated.  Function  subprograms 
were  used  to  perform  summing  operations  and  computation 
of  the  variance  of  the  a parameter  due  to  its  complexity. 
Table  £.1  lists  the  variable  names  used  in  the  main 
program  with  their  associated  meanings.  Figure  £.3  is  a 
simplified  flowchart  of  the  main  program. 
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Table  £.1 

VARIABLES  USED  IV  LXKELX2 


AVS 

Character  variable  uaed  to  determine  If 
full  explanation  of  program  usage  is 
desired. 

FILEKAME 

Character  variable  used  to  contain  the 
name  of  the  file  oontaining  data  to  be 
analysed. 

JR 

Counter  used  in  summing  subprograms. 

J 

Counter  uaed  in  do-loop. 

T(J) 

Vector  used  to  store  failure  data  points. 

XM 

The  value  for  the  Weibull  parameter  "m" 
for  which  this  program  was  developed. 

BIGT 

The  truncation  point  specified. 

N 

Humber  of  elements  in  original  population 

X 

Individual  data  points * as  read  from  file 

I 

Counter  used  to  determine  number  of  data 
points  read  in. 

ALT 

Variable  used  to  store  the  computed  value 
of  the  left  side  of  Equation  (E.3). 

ART 

Variable  used  to  store  the  computed  value 
of  the  right  side  of  Equation  (E.3). 

XK 

The  value  for  the  Weibull  parameter  MK" 
for  which  this  program  was  developed. 

ZZ 

Dummy  variable  used  in  reading  the  line 
numbers  on  the  data  file  (once  read  these 
numbers  are  discarded). 

XLOW 

Variable  used  to  store  the  value  of  XK 
for  interval  bisection  computations. 
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XHIGH  Variable  used  to  store  the  value  of  XM 

for  interval  bisection  computations. 

XDIF  Difference  between  XHIGH  and  XLOW. 

XKVAR  Variance  of  Weibull  rtKw  parameter  esti- 

mate computed  using  MLE’s. 
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below  ALT  when  the 
program  is  -i  uses 
interval  bisect  as  in 

1. 


Figure  E.3 

Giaplif'ed  Flowchart  of  Main  Program 
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mums  fork  OF  Mas  algorithm 


APP9DIX  9 

fobthai  worn  or  use  Auxmvm 

The  algorithm  developed  in  Chapter  IV  was 
pro  grassed  in  PORTKAV  and  is  shown  in  Figure  ?.?. 

basically*  the  program  consists  of  a Main  Program  and 

\ ■ 

five  function  subprograms.  The  main  program  initialises 
failure  distribution  parameters  and  cost  inputs*  contains 
two  do-loops  which  increment  the  value  of  n and  the  age 
of  the  core  module*  manipulates  probabilities  and  expec- 
tations returned  by  the  subprograms*  and  provides  data 
output.  Tho  function  subprograms  perform  the  detailed 
calculations  of  probabilities  and  eaqpected  values. 

Table  F.i  lists  the  variable  names  used  in  the 
main  program  with  their  associated  meanings.  Throughout 
the  program,  MODI  refers  to  the  module  for  which  the 
vp  .-al  opportunistic  replacement  policy  is  being  deter- 
s.  ..  < d.  M0D2  refers  to  the  core  module.  The  program  could 
be  used  for  any  two  component  system  where  it  is  desired 
to  find  the  optimal  opportunistic  replacement  policy  for 
one  of  the  components  given  the  maximum  operating  times 
and  failure  distribution  parameters.  For  this  reason* 
the  symbols  MODI  and  MOSS  were  used  rather  than*  for 
instance,  fan  and  core. 
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Numerical  solution  for  the  optimal  vmlue  of  n ie 
performed  in  the  main  program.  Figure  P.1  ie  a simpli- 
fied flow  chert  of  the  main  program.  The  value  of  MOD2AGE 
ic  initially  set  at  0 at  line  400.  The  value  of  TSTAB  is 
initially  set  to  0 at  line  600.  Lines  610  through  640 
determine  the  difference*  if  any*  between  M0D2AGE  and 
TSTAB,  calculate  the  maximum  remaining  hours  on  MODI, 
and  determine  if  the  value  of  N represents  MODI  MOT  or 
M0D2M0T.  Lines  660  through  840  make  the  necessary  calls 
to  the  function  subprograms  to  determine  values  of  CYCLE, 
TCQST  and  HOUKCOST.  These  three  values  are  printed  out 
along  with  M0D2AGE  and  TSTAB.  The  program  then  loops 
back,  increments  the  value  of  TSTAB  by  MOD1MOT/25  hours 
and  determines  CYCLE,  TCOST  and  HOUBCOST  once  more.  TSTAR 
is  incremented  in  MODI MOT/25  hour  steps  until  MODI MOT  is 
reached.  Once  M0D1M0T  is  reached,  the  program  loops  back 
to  line  440,  increments  M0D2AGE  by  the  value  of  SCAN, 
sets  the  value  of  TSTAB  to  zero,  and  increments  TSTAB 
in  MODiMOT/25  hour  steps  once  more.  This  pattern  ol‘ 
increment ing  TSTAB  from  0 hours  to  MOT  hours  and  then 
incrementing  the  value  of  M0D2AGE  is  continued  until 
M0D2AGE  equals  25  times  the  value  of  SCAN.  If  SCAN  is  set 
at  M0D2M0T/25*  the  final  value  of  M0D2AGE  will  equal 
M0D2M0T. 

As  currently  written,  the  program  prints  one  line 
of  output  for  each  value  of  M0D2AGE  and  TSTAB.  By  deleting 
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line  9/0  and  inserting  the  following  statement  between  lines 
930  and  9^0,  the  program  will  print  the  optimal  value  of 
TSTAR  (and  the  associated  values  of  CYCLE*  TCOST  and  HOUR- 
COST)  for  each  value  of  M0D2AGE: 

PRINT  5*  M0D2AGE ,XSTAR *H0URC0ST  * CYCLE  * TCOST 

Function  subprogram  SINGMEAN  performs  the  integral 

30,01  • -T61**  t . *M0D1  <t> « (,.d 

numerically  using  Simpson* s rule  (23:14).  Numerical 
evaluation  was  necessary  since  (F.1)  has  a closed  form  for 
the  Weibull  PDF  only  when  the  limits  of  integration  are  0 
and  infinity.  In  order  to  insure  that  SINGMEAN  functioned 
properly,  a special  "test'*  Main  Program  was  designed  with 
Just  sufficient  statements  to  pass  arguments  to  SINGMEAN 
and  print  results.  The  Welbull  PDF  was  sketched  for 
m * .71019  and  k • ,000163  which  were  the  maximum  likelihood 
estimates  of  the  parameters  of  the  fan  drive  turbine  module. 
The  value  of  the  PDF  is  approximately  zero  for  values  of  t 
beyond  300  hours.  Mean  time  to  failure  was  analytically 
determined  as  (24:221) 

» • (a-pr)  j 

The  analytical  mean  of  a Weibull  distribution  with 
m ■ .71019  and  k ■ .000163  is  200.172  hours.  A value  of 
t « 6000  hours  which  is  well  beyond  300  hours  was  passed 
to  SINGMEAN.  The  absolute  error  of  the  value  of  XTHJ1 
returned  by  SINGMEAN  is  sensitive  to  the  number  of  segments 
used  in  numerical  approximation  of  the  integral.  With  200 
segments*  error  was  over  six  percent.  By  increasing  the 


IWSWW sir* 1 m - . - 


number  of  segments  to  2000 « the  am  returned  fey 
mi  201.683  which  is  Hjwiflwtily  .6  percent  error.  lbs 
relatively  large  master  of  segments  used  in  BTKitltH  is 
the  primary  reason  the  algorithm  requires  1.2  hoars  of 
processing  tine.  Further  increasing  the  amber  of  seg- 
ments in  SHGMSAV  would  require  even  longer  processing 
time.  Tor  this  reason*  .8  percent  error  was  accepted. 

The  complete  algorithm  was  run  once  with  3000  segments 
to  determine  if  the  optimal  value  of  TSTAfi  was  sensitive 
to  residual  error  in  SZffGfCEAV.  The  optimal  value  of  T8TAR 
wa a the  same  aa  that  found  using  2000  segments  in  SZTOHSHih. 
HOURCOST  was  approximately  .07  percent  higher  using  3000 
segments.  This  difference  is  felt  to  be  insignificant. 

Funetion  subprogram  QWOD  calculates  the  value  of 


. FCt)  - F(t,) 


(P.?> 


for  the  Weibull  PDF.  No  approximations  are  used  in  this 
function  subprogram.  The  subprogram  was  checked  individually* 
however*  in  the  same  manner  as  SINGMEAN  and  returned  the 
analytically  correct  probability. 

Function  subprograms  QFCEAR*  QFBARC  and  QQFC  are 
essentially  similar.  The  mechanics  of  the  three  function 
subprograms  are  illustrated  here  by  reference  to  function 
QFCBAR.  In  Chapter  IV*  the  probability  of  the  fan  module 
requiring  replacement  before  the  core  module  in  the 
interval  (n*  N)  given  fan  module  survival  until  n and  core 


i.£o 


i m 1 1TT 


module  survival  until  n • 4t  was 


M 

Q(fc)  - £ *!<*«)  (*.*) 

i-1 


when  n<  t^H  - It, 
where 


- [l  - «PAN(ti  | !>)]• 

[l  - QOORE(ti+1  | n)]. 


(F.5) 


|ti) 

and  when  N - I t <ttSS  q-^fc)  - 0 if  N - COROIOT  but  if 


N - FAHMOT  then 


f2D  - [l  - QFAK(»  - I t | n)]e  (F.6) 

[l  - QCORE(N  | n)J. 

Function  QFCBAR  performs  this  summation.  Three  values  are 
returned  by  the  subprogram.  In  addition  to  the  value  of 
(F.3)  the  subprogram  also  returns 

H-1  (t  * t ) 

QFCBARL  - £ ^ 1 * titV  q.(fc)  (1*’.7) 

iTi  5 


and 


M-1  <t,  * t.  . i 

imn  - £ — q. as)  U'.e; 

i.i 

QFCBARL  is  used  to  pass  the  value  of  XA  to  the  main  program. 
XHtT>  is  used  to  pass  the  contribution  of  Q(fc)  to  the  cal- 
culation of 
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£ C*i  ♦ f<,.(fc)  + q1(Fo>+qi(fo)|  (F.9) 

i.i  2 L J 

required  by  Equation  (4.17)  in  Chapter  IV.  Function  sub- 
programs QftFC  and  QFBARC  differ  from  function  QFCBAR  only 
in  that  they  perform  calculations  involving  qi(fc)  and  qi(Fc) 

respectively. 

Function  subprograms  QFCBAR,  QFBARC  and  QQFC  were 
checked  toy  writing  a short  main  program  with  just  sufficient 
statements  to  pass  the  necessary  arguments  and  print  results. 
The  value  of 


M 


1 - £ q^fc)  + qi(7c)  + qi(fc) 
i*1 


(F.10) 


over  the  interval  (0,  t)  is  equivalent  to 
['  - »*»!<*>]  * ['  - 

Thus  comparing  the  value  of  (F.10>  with  the  analytic  two 
component  system  reliability  for  some  time  t provides  a 
check  on  the  function  subprograms.  This  check  was  per- 
formed for  a value  of  XM1  * .66,  XK1  - .000141,  XM2  = -.1316 
and  XK2  * .002051.  The  analytic  reliability  is  .334.  With 


200  segments  In  each  of  the  function  subprograms,  the  value 
of  (F.10)  was  .336  which  is  .6  percent  error. 

A second  check  of  subprograms  QFCBAR,  QFBARC  and 
QftFC  is  possible  by  comparing  the  mean  time  at  failure 
of  a two  component  system  returned  by  the  subprograms 
with  the  analytically  determined  mean  time  at  failure. 
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For  a system  of  Weibull  components*  mean  time  at  failure 
may  be  calculated  an  .lytically  over  the  interval  (0*  ••  ) 

A 

when  the  component  m parameters  are  equal  (23:221 ) For 
m = .66004  and  k values  of  .000141  and  .002051*  analytic 
mean  time  at  failure  is  48.38  hours.  The  m and  k values 
chosen  are  in  the  range  relevant  to  distributions  studied 
in  this  thesis.  The  mean  determined  by  QFCBAW , QFBAHC 
and  QQFC  was  48.48  hours  which  is  an  error  of  approximately 
.2  percent.  Functions  QFCBAR*  QFBARC  and  QQFC*  like 
SINGMEAN*  require  an  upper  bound  for  the  interval  over 
which  approximation  is  to  be  accomplished.  For  the  m or 
k values  chosen*  the  PDF's  are  approximately  0 beyond 
t * 500  hours.  An  upper  bound  of  2000  was  used  to  check 
the  subprograms. 


'The  formula  for  calculating  system  mean  time  to 
failure  is 


■ ^tT)*Pm+1 }(  a + ,?  ) j 


i/(m+i  j 


Table  P.1 
PROGRAM  SYMBOLOGY 
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N 

MOD2AGE 
MODI PACK 

MODI DEP 
MODI REP 


MODI SHI? 


MODI MOT 


MGD2M07 


;/v- 


XK' 

m 

XK? 

SCAIJ 


Maxinun  possible  age  at  replacement  of 
module  1.  Refer  to  Equation  (4.21 ) in 
Chapter  IV. 

Age  of  module  2. 

Dollar  cost  to  pack  and  unpack  module  1 at 
the  field  level. 

Dollar  cost  to  depot  overhaul  module  1. 

Dollar  cost  to  remove  and  replace  module  1. 

This  cost  includes  the  cost  to  remove  the 
engine  from  the  aircraft,  transport  the 
engine  to  the  repair  shop*  perform  the  in 
shop  module  removal  and  replacement,  trans- 
port the  engine  from  the  shop  to  the  aircraft;, 
reinstall  and  trim  check  the  engine. 

Cost,  in  dollars,  to  transport  module  1 between 
the  field  repair  shop  and  the  depot.  Shipping 
cost  is  round  trip. 

Module  1 maximum  operating  time. 

Module  2 maximum  operating  time. 

Weibull  m parameter  for  module  ;. 

W'  i.buli  X parameter  for  module  1. 

Weibuli.  n parameter  for  module  ; . 

Weibull  k parameter  for  module  2. 

Amount  oy  wnich  module  2 age  is  incremented 
in  d j-loop. 

Di-  erence  between  module  ' MOT  and  module 
2 MOT.  Refer  to  Equation  'h.20)  in  Chapter 


DELTAM 


■ww*  '.iii'  i#b  .wjjCTTBagwra  i *■ 1 y^ffapgyy,w^vnv'? 
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TSTAR 

TCOST 

CYCLE 

HOURCOST 

XLOW 

XSTAR 
I.  J 
UL 

DELTAT 

NPIX 

XI 


X? 


X3 


X4 

X5 

X6 

X7 


Breakpoint  between  replace  at  failure  region 
and  opportunistic  replacement  region.  TSTAR 
is  equivalent  to  n as  used  in  Chapter  IV. 

Total  cost  per  cycle  in  dollars. 

Cycle  length  in  hours. 

Cost  per  hour  in  dollars  per  hour. 

Variable  used  to  store  the  lowest  value  of 
hourcost . 

The  value  of  TSTAR  associated  with  XLOW. 

Counters  used  in  do-loops. 

Maximum  hours  remaining  on  module  1 before 
replacement . 

Difference  between  module  1 age  and  module 
2 age.  Refer  to  Equation  (4.19)  in  Chapter 
IV. 

Switch  variable  used  to  pass  information  on 
which  module  will  reach  its  MOT  first. 

Probability  of  module  1 requiring  replacement 
before  module  2 in  the  interval  n<ts;N 
given  survival  of  both  until  TSTAR  and  TSTAR  - 
DELTAT  respectively. 

Probability  of  module  2 requiring  replacement 
before  module  1 in  the  interval  TSTAR ■< t sN 
given  survival  of  both  until  TSTAR  and 
TSTAR- DELTAT  respectively. 

jProbability  of  module  i and  module  2 requiring 
simultaneous  replacement  in  the  interval 
TSTAR  <t  sN  given  survival  of  both  until 
TSTAR  and  TSTAR-  DELTAT  respectively. 

Defined  analogously  to  XI  over  the  interval 
TSTAR  «t  <N. 

Defined  analogously  to  X2  over  the  interval 
TSTAR «*t  «N. 

Defined  analogously  to  X3  over  the  interval 
TSTAR  < t «R. 

Probability  of  module  1 failure  in  the  interval 
O^taS  TSTAR. 


Be liability  of  module  1 in  the  interval 


Probability  that  module  1 requires  replace- 
ment before  module  ? in  the  interval 
TSTAR  <tjS N given  module  ? survival  until 
TSTAR- DELTAT. 

Probability  that;  module  requires  ^place- 
ment before  module  1 in  the  interval 
TSTAR  ^tjfiN  given  module  ? survival  until 
TSTAR- DELTAT . 

Probability  that  module  1 and  module  ? require 
simultaneous  replacement  in  the  interval 
TSTAR  «t  *N  given  module  2 survival  until 
TSTAR-DELTAT. 

Defined  analogously  to  XII  over  the  interval 
TSTAR  <*srt  <N. 

Defined  analogously  to  X22  over  the  interval 
TSTAR  <e  t <tN. 

Defined  anal  ^ously  to  X33  over  the  interval 
TSTAR  «t«N. 


The  contribution  to  mean  time  to  removal  of 
module  1 by  the  interval  OjftaTSTAR. 

The  contribution  to  mean  time  to  removal  of 
module  1 by  the  interval  TSTAR  <t<N. 
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Figure  t.A 

Sioplified  Flowchart  of  Main  Prograa  H3SWL-$frU 
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AFFODH  G 


F«X»AH  OUTPUT 

This  appendix  contains  output  data  from  the 
FOBTSAH  program  listed  in  Appendix  F.  Figure  G.l  is  a 
8 ample  of  the  output  provided  by  the  program.  Table  G.l 
summarises  the  results  from  sensitivity  analysis  and  is 
presented  in  the  following  formats 

1.  Column  1 indicates  the  input  variable  under  study. 

2.  Column  ? contains  the  value  of  the  input  variable 
under  study. 

3.  The  value  in  column  3 is  the  optimal  value  of  n 
(the  breakpoint  between  the  replace  at  failure  region  add 
the  opportunistic  replacement  region)  for  the  input  para- 
meters used.  Precision  is  +,  10  hours. 

4.  The  value  in  column  4 is  the  expected  cost  per 
hour  under  the  optimal  opportunistic  replacement  policy 
for  the  input  parameters  used. 

5.  The  value  in  column  5 is  the  expected  cost  per 
hour  under  the  corresponding  replace  at  failure  policy. 

6.  The  value  in  column  6 is  the  difference  in 
expected  cost  between  the  replace  at  failure  policy  and 
the  opportunistic  replacement  policy. 

7.  The  value  in  column  7 is  the  expected  cycle  length 
under  the  optimal  opportunistic  policy  for  the  input 
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parameter  used. 

8.  The  value  in  column  8 is  the  expected  cycle  length 
under  the  corresponding  replace  at  failure  policy.  The 
baseline  values  of  the  program  input  variables  were: 


XM1 

.71019 

XK1 

.000165 

xmi 

-.01516 

XK2 

.002051 

MODI REP 

$862.52 

MODI DEP 

$61 7^ • 00 

MODI PACK 

$50.10 

MODI SKIP 

$315.84 

MODI MOT 

250.0  hours 

M0D2M0T 

250.0  hours 

SC  All 

10.0  hours 

Table 


1?1 


M0D2AGE 


10.0000 
10.  0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

•10.0000 

10.0000 
"0.0000 

10.0000 
"0.0000 

10.0000 

10.0000 

10.0000 

10.0000 


TSTAfi 


0. 

10.0000 

70.0000 

50.0000 

40.0000 

50.0000 

60.0000 

70.0000 

80.0000 

90.0000 

100.0000 

110.0000 

120.0000 

150.0000 

140.0000 

150.0000 

160.0000 

170.0000 

180.0000 

190.0000 

200.0000 

210.0000 

220.0000 

230.0000 

240.0000 

250.0000 


*»  ‘IHCOST 


49.**'  -*i 

49.550 

48.791 

48.097' 

47.46 

46.88s> 

46.3730 

45.9091 

45.4929 

45.  "1206 

44.7889 

44.4948 

44.2355 

44.0084 

43. 8111 

43.6416 

45.4978 

43.3779 

43. 2803 

43.2034 

43.1458 

43.1004 

45.0838 

43.0771 

45. 0853 

43.1077 


CYCLE 


139.1934 
143.8107 
146. 5348 
148.7167 
150. 9506 
153.0554 
155.0147 

156.8368 
158.5243 
160.0803 
161.5084 
162.8124 
163.9964 
165.0648 
166.0220 
166.8725 
167.6212 
168. 2725 
168.8312 
169.3020 
169.6892 
109.9974 
170.2508 
170.3935 
170.4895 
170.5224 


Figure  G.i 

Sample  Program  Output ' 


A t 

baa*; lino  hazard  rate  and  cost  inputs 
dr Lve  turbine;  Oare  Age  » 10  hours. 


TCOST 


6938.6748 
7125.9541 
7139.86 27 
7152.1002 

7104.4039 

7176.5.373 

7188.4919 

7700.2351 
771 1 . 721 1 
7227.9142 
7733.7845 
7244.3069 
7254.4614 

7264.2319 
7273.6074 
7282.5801 
7291.1464 
7299.3057 
7307.0007 
7514.4  171 
737--.5820 
7327.9084 
7354.1871 
7 340.0540 
7345. 5899 
7350.6700 
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CONDITIONAL  PROBABILITIES 


APPEBD1X  H 

COVDITIOKAL  PROBABILITIES 

Chapter  IV  made  extensive  use  of  conditional 
probabilities  of  the  fora  q^(fc)»  q^(?c)  and  q^(fd’). 

The  reader  will  recall  * for  ins  tance*  that 

qi(fc)  - - QPAN(t.  | n)  ja  (H.1 ) 

Tl  - QCORE(ti+1J| 

QPAN(ti+1  | t.)  • 

*;oi<K(ti+1 1 t.  ) 

ir.  the  internal  n<  t^<  N - I t , 


where 

q.(fc)  ■ Probability  that  an  engine  requires 
1 field  level  replacement  of  the  fan 

module  and  core  module  simultaneously 
in  the  i-th  interval  between  n and  N 
given  survival  of  both  until  n*  and 

^FAUCt^  It,  ) a 
“•Fan  J dFan 

Probability  of  fan  failure  between  f . 

U 

and  t.  given  survival  until  t .. 

•.*OORI*.C  t.  It,  ) ~ 

“core  I “core 

Probability  of  core  failure  between  t, 

o 

and  given  survival  until  t ^ . 

Let 

fpan(t;IFan  module  failure  density  function. 

f (t)*Core  module  failure  density  function, 
core 
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which  reduces  to 


ir.  ilarly , 


n)j. 


^CORh(t 
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and  Equation  (H.l)  becomes 

q^fc)  . (H.3) 


Equation  (H.3)  represents  the  probability  ef  FAN/COltE 
system  failure  in  the  interval  (t^*  t^+1)  given  survival 

until  n.  The  devisor  in  Equation  (H.3)  insures  that  the 
probability  of  failure  over  the  interval  (n,®*»)  is 
equal  to  1 . A geometric  interpretation  of  the  conditional 
probabilities  in  Equation  (H.3)  is  shpwn  in  Figure  H.l. 

The  PDF  has  a spike  at  N where  mandatory  removal  occurs. 
q^Tc)  and  q^(fc)  are  developed  analogously  with  the 

exception  that  one  module  survives  through  the  interval 
^ t i , ti+1)  while  the  other  fails  in  it. 
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